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ABSTRACT: Cosmic microwave background (CMB) temperature anisotropies have and will 
continue to revolutionize our understanding of cosmology. The recent discovery of the previously 
predicted acoustic peaks in the power spectrum has established a working cosmological model: 
a critical density universe consisting of mainly dark matter and dark energy, which formed its 
structure through gravitational instability from quantum fluctuations during an inflationary 
epoch. Future observations should test this model and measure its key cosmological parameters 
with unprecedented precision. The phenomenology and cosmological implications of the acoustic 
peaks are developed in detail. Beyond the peaks, the yet to be detected secondary anisotropies 
and polarization present opportunities to study the physics of inflation and the dark energy. The 
analysis techniques devised to extract cosmological information from voluminous CMB data sets 
are outlined, given their increasing importance in experimental cosmology as a whole. 

1 INTRODUCTION 

The field of cosmic microwave background (CMB) anisotropies lias dramatically 
advanced over the last decade (c.f. White et al 1994), especially on its obser- 
vational front. The observations have turned some of our boldest speculations 
about our Universe into a working cosmological model: namely, that the Universe 
is spatially fiat, consists mainly of dark matter and dark energy, with the small 
amount of ordinary matter necessary to explain the light element abundances, 
and all the rich structure in it formed through gravitational instability from quan- 
tum mechanical fluctuations when the Universe was a fraction of a second old. 
Observations over the coming decade should pin down certain key cosmological 
parameters with unprecedented accuracy (Knox 1995; Jungnian et al 1996; Bond 
et al 1997; Zaldarriaga et al 1997; Eisenstein et al 1999). These determinations 
will have profound implications for astrophysics, as well as other disciplines. Par- 
ticle physicists, for example, will be able to study neutrino masses, theories of 
inflation impossible to test at accelerators, and the mysterious dark energy or 



cosmological constant. 

For the twenty eight years between the discovery of the CMB (Penzias & Wil- 
son 1965) and the COBE DMR detection of 10"'"' fluctuations in its temperature 
field across the sky (Snioot et al 1992), observers searched for these anisotropies 
but found none except the dipole induced by our own motion (Smoot et al 1977). 
They learned the hard way that the CMB is remarkably uniform. This is in stark 
contrast to the matter in the Universe, organized in very non-linear structures 
like galaxies and clusters. The disparity between the smooth photon distribution 
and the clumpy matter distribution is due to radiation pressure. Matter inho- 
mogeneities grow due to gravitational instability, but pressure prevents the same 
process from occuring in the photons. Thus, even though both inhomogeneities 
in the matter in the Universe and anisotropies in the CMB apparently originated 
from the same source, these appear very different today. 

Since the photon distribution is very uniform, perturbations are small, and 
linear response theory applies. This is perhaps the most important fact about 
CMB anisotropies. Since they are linear, predictions can be made as precisely 
as their sources are specified. If the sources of the anisotropies are also linear 
fiuctuations, anisotropy formation falls in the domain of linear perturbation the- 
ory. There are then essentially no phenomenological parameters that need to be 
introduced to account for non-linearities or gas dynamics or any other of a host 
of astrophysical processes that typically afflict cosmological observations. 

CMB anisotropies in the working cosmological model, which we briefly review 
in §2, fall almost entirely under linear perturbation theory. The most important 
observables of the CMB are the power spectra of the temperature and polariza- 
tion maps. Theory predicts, and now observations confirm, that the temperature 
power spectrum has a series of peaks and troughs. In §3, we discuss the origin 
of these acoustic peaks and their cosmological uses. Although they are the most 
prominent features in the spectrum, and are the focus of the current generation 
of experiments, future observations will turn to even finer details, potentially re- 
vealing the physics at the two opposite ends of time. Some of these are discussed 
in §4. Finally, the past few years have witnessed important new advances, intro- 
duced in §5, from a growing body of CMB data analysts on how best to extract 
the information contained in CMB data. Some of the fruits of this labor have 
already spread to other fields of astronomy. 

2 OBSERVABLES 

2.1 Standard Cosmological Paradigm 

While a review of the standard cosmological paradigm is not our intention (see 
Narkilar & Padmanabhan 2001 for a critical appraisal), we briefiy introduce the 
observables necessary to parameterize it. 
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The expansion of the Universe is described by the scale factor a{t), set to unity 
today, and by the current expansion rate, the Hubble constant Hq = lOO/i km 
sec~^ Mpc"'^, with h ~ 0.7 (Preedman et al 2001). The Universe is flat (no spatial 
curvature) if the total density is equal to the critical density, pc = 1.88/i^ x 10~^^g 
cm""'*; it is open (negative curvature) if the density is less than this and closed 
(positive curvature) if greater. The mean densities of different components of the 
Universe control a{t) and are typically expressed today in units of the critical 
density Oj, with an evolution with a specified by equations of state Wj = Pi/pi, 
where pi is the pressure of the ith component. Density fluctuations are determined 
by these parameters through the gravitational instability of an initial spectrum 
of fluctuations. 

The working cosmological model contains photons, neutrinos, baryons, cold 
dark matter and dark energy with densities proscribed within a relatively tight 
range. For the radiation, Qj. = 4.17 x 10~^h~^ {wj. = 1/3). The photon contri- 
bution to the radiation is determined to high precision by the measured CMB 
temperature, T = 2.728 ±0.004K (Fixsen et al 1996). The neutrino contribution 
follows from the assumption of 3 neutrino species, a standard thermal history, 
and a negligible mass <^ leV. Massive neutrinos have an equation of state 
= 1/3 — » as the particles become non-relativistic. For ~ leV this occurs 
at a 10~^ and can leave a small but potentially measurable effect on the CMB 
anisotropies (Ma & Bertschinger 1995; Dodelson et al 1996). 

For the ordinary matter or baryons, flf, « 0.02h~^ (wt ~ 0) with statistical 
uncertainties at about the ten percent level determined through studies of the 
light element abundances (for reviews, see Bocsgaard & Stcigman 1985; Schramm 
& Turner 1998; Tytler et al 2000). This value is in strikingly good agreement 
with that implied by the CMB anisotropies themselves as wc shall see. There 
is very strong evidence that there is also substantial non-baryonic dark matter. 
This dark matter must be close to cold {wm = 0) for the gravitational instability 
paradigm to work (Peebles 1982) and when added to the baryons gives a total 
in non-relativistic matter of n,„ — 1/3. Since the Universe appears to be flat, 
the total Qtot must be equal to one. Thus, there is a missing component to 
the inventory, dubbed dark energy, with Ua ~ 2/3. The cosmological constant 
(wA = — 1) is only one of several possible candidates but we will generally assume 
this form unless otherwise specified. Measurements of an accelerated expansion 
from distant supcrnovac (Riess ot al 1998; Perlmutter et al 1999) provide entirely 
independent evidence for dark energy in this amount. 

The initial spectrum of density perturbations is assumed to be a power law 
with a power law index or tilt of n ~ 1 corresponding to a scale-invariant spec- 
trum. Likewise the initial spectrum of gravitational waves is assumed to be scale- 
invariant, with an amplitude parameterized by the energy scale of inflation Ei, 
and constrained to be small compared with the initial density spectrum. Finally 
the formation of structure will eventually reionize the Universe at some redshift 
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7 < 2ri ;S 20. 

Many of the features of the anisotropies will be produced even if these parame- 
ters fall outside the expected range or even if the standard paradigm is incorrect. 
Where appropriate, we will try to point these out. 

2.2 CMB Tem.perature Field 

The basic observable of the CMB is its intensity as a function of frequency and 
direction on the sky n. Since the CMB spectrum is an extremely good blackbody 

(Fixsen et al 1996) with a nearly constant temperature across the sky T, we 
generally describe this observable in terms of a temperature fluctuation B(n) = 
AT/T. 

If these fluctuations are Gaussian, then the multipole moments of the temper- 
ature field 

Qem = ldhY;^{n)Q{n) (1) 
are fully characterized by their power spectrum 

{®em®e'm') = See'Smm'Ce , (2) 

whose values as a function of £ are independent in a given realization. For this 
reason predictions and analyses are typically performed in harmonic space. On 
small sections of the sky where its curvature can be neglected, the spherical 
harmonic analysis becomes ordinary Fourier analysis in two dimensions. In this 
limit £ becomes the Fourier wavcnumber. Since the angular wavelength $ = 
27r/^, large multipole moments corresponds to small angular scales with £ ~ lO'^ 
representing degree scale separations. Likewise, since in this limit the variance 
of the field is / d^£Q/(27r)^, the power spectrum is usually displayed as 

^2^^e{^CeT\ (3) 

the power per logarithmic interval in wavenumber for £ 2> 1. 

Plate 1 (top) shows observations of At' along with the prediction of the work- 
ing cosmological model, complete with the acoustic peaks mentioned in §1 and 
discussed extensively in §3. While COBE first detected anisotropy on the largest 
scales (inset), observations in the last decade have pushed the frontier to smaller 
and smaller scales (left to right in the figure). The MAP satellite, launched in 
June 2001, will go out to £ ^ 1000, while the European satellite. Planck, sched- 
uled for launch in 2007, will go a factor of two higher (see Plato 1 bottom). 

The power spectra shown in Plate 1 all begin at £ = 2 and exhibit large 
errors at low multipoles. The reason is that the predicted power spectrum is the 
average power in the multipole moment £ an observer would see in an ensemble 
of universes. However a real observer is limited to one Universe and one sky with 
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its one set of G^m's, 2£ + 1 numbers for each £. This is particularly problematic 
for the monopole and dipole {I = 0,1). If the monopole were larger in our 
vicinity than its average value, we would have no way of knowing it. Likewise 
for the dipole, we have no way of distinguishing a cosmological dipole from our 
own peculiar motion with respect to the CMB rest frame. Nonetheless, the 
monopole and dipole - which we will often call simply O and v-y - are of the utmost 
significance in the early Universe. It is precisely the spatial and temporal variation 
of these quantities, especially the monopole, which determines the pattern of 
anisotropies we observe today. A distant observer sees spatial variations in the 
local temperature or monopole, at a distance given by the lookback time, as a 
fine-scale angular anisotropy. Similarly, local dipoles appear as a Doppler shifted 
temperature which is viewed analogously. In the jargon of the field, this simple 
projection is referred to as the freestreaming of power from the monopole and 
dipole to higher multipole moments. 

How accurately can the spectra ultimately be measured? As alluded to above, 
the fundamental limitation is set by "cosmic variance" the fact that there are 
only 2£ + 1 m-samples of the power in each multipole moment. This leads to an 
inevitable error of 

AC. = ^C.. (4) 

Allowing for further averaging over t in bands of A£ ~ £, we see that the precision 
in the power spectrum determination scales as £~^, i.e. ~ 1% at £ = 100 and 
~ 0.1% a,t I = 1000. It is the combination of precision predictions and prospects 
for precision measurements that gives CMB anisotropies their unique stature. 

There are two general caveats to these scalings. The first is that any source of 
noise, instrumental or astrophysical, increases the errors. If the noise is also Gaus- 
sian and has a known power spectrum, one simply replaces the power spectrum 
on the rhs of Equation (4) with the sum of the signal and noise power spectra 
(Knox 1995). This is the reason that the errors for the Planck satellite increase 
near its resolution scale in Plate 1 (bottom). Because astrophysical foregrounds 
are typically non-Gaussian it is usually also necessary to remove heavily contam- 
inated regions, e.g. the galaxy. If the fraction of sky covered is /sky, then the 

— 1/2 

errors increase by a factor of /^.j^^ and the resulting variance is usually dubbed 
"sample variance" (Scott et al 1994). An /s^y = 0.65 was chosen for the Planck 
satellite. 

2.3 CMB Polarization Field 

While no polarization has yet been detected, general considerations of Thomson 
scattering suggest that up to 10% of the anisotropies at a given scale are polar- 
ized. Experimenters are currently hot on the trail, with upper limits approaching 
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Plate 1: Top: temperature anisotropy data with boxes representing l-a errors and approximate 
^-bandwidth. Bottom: temperature and polarization spectra for Otot — 1, ^^A — 2/3, Qbh^ — 
0.02, n,„/i^ = 0.16, n = 1, Zri = 7, Ei ^ 2.2 X lO^** GeV. Dashed lines represent negative cross 
correlation and boxes represent the statistical errors of the Planck satellite. 
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the expected level (Hedman et al 2001; Keating et al 2001). Thus, we expect 
polarization to be an extremely exciting field of study in the coming decade. 

The polarization field can be analyzed in a way very similar to the temper- 
ature field, save for one complication. In addition to its strength, polarization 
also has an orientation, depending on relative strength of two linear polarization 
states. While classical literature has tended to describe polarization locally in 
terms of the Stokes ijaramcters Q and U^, recently cosmologists (Scljak 1997; 
Kamionkowski et al 1997; Zaldarriaga & Seljak 1997) have found that the scalar 
E and pseudo-scalar B, linear but non-local combinations of Q and U, provide 
a more useful description. Postponing the precise definition of E and B until 
§3.7, we can, in complete analogy with Equation (1), decompose each of them 
in terms of multipole moments, and then, following Equation (2), consider the 
power spectra, 

i^lm^e'm') = ^U'^mm'Cf^ ■ (5) 

Parity invariance demands that the cross correlation between the pseudoscalar B 
and the scalars or _E vanishes. 

The polarization spectra shown in Plate 1 [bottom, plotted in fiK following 
Equation (3)] have several notable features. First, the amplitude of the EE spec- 
trum is indeed down from the temperature spectrum by a factor of ten. Second, 
the oscillatory structure of the EE spectrum is very similar to the temperature 
oscillations, only they arc apparently out of phase but correlated with each other. 
Both of these features are a direct result of the simple physics of acoustic oscil- 
lations as will be shown in §3. The final feature of the polarization spectra is 
the comparative smallness of the BB signal. Indeed, density perturbations do 
not produce B modes to first order. A detection of substantial B polarization, 
therefore, would be momentous. While E polarization effectively doubles our cos- 
mological information, supplementing that contained in Ce, B detection would 
push us qualitatively forward into new areas of physics. 



3 ACOUSTIC PEAKS 

When the temperature of the Universe was ~ 3000K at a redshift « 10^, 
electrons and protons combined to form neutral hydrogen, an event usually known 

as recombination (Peebles 1968; Zel'dovich et al 1969; see Seager et al 2000 
for recent refinements). Before this epoch, free electrons acted as glue between 
the photons and the baryons through Thomson and Coulomb scattering, so the 

^There is also the possibility in general of circular polarization, described by Stokes parameter 
V, but this is absent in cosmological settings. 
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cosmological plasma was a tightly coupled photon-baryon fluid (Peebles & Yu 
1970). The spectrum depicted in Plate 1 can be explained almost completely by 
analyzing the behavior of this pre-recombination fluid. 

In §3.1, we start from the two basic equations of fluid mechanics and derive 
the salient characteristics of the anisotropy spectrum: the existence of peaks 
and troughs; the spacing between adjacent peaks; and the location of the first 
peak. These properties depend in decreasing order of importance on the initial 
conditions, the energy contents of the Universe before recombination and those 
after recombination. Ironically, the observational milestones have been reached in 
almost the opposite order. Throughout the 1990"s constraints on the location of 
the first peak steadily improved culminating with precise determinations from the 
TOCO (MUler et al 1999), Boomerang, (de Bernardis et al 2000) and Maxima- 1 
(Hanany et al 2000) experiments (see Plate 1 top). In the working cosmological 
model it shows up right where it should be if the present energy density of the 
Universe is equal to the critical density, i.e. if the Universe is flat. The skeptic 
should note that the working cosmological model assumes a particular form for 
the initial conditions and energy contents of the Universe before recombination 
which we shall see have only recently been tested directly (with an as yet much 
lower level of statistical confidence) with the higher peaks. 

In §3.2 we introduce the initial conditions that apparently arc the source of all 
dumpiness in the Universe. In the context of ab initio models, the term "initial 
conditions" refers to the physical mechanism that generates the primordial small 
perturbations. In the working cosmological model, this mechanism is inflation 
and it sets the initial phase of the oscillations to be the same across all Fourier 
modes. Remarkably, from this one fact alone comes the prediction that there 
will be peaks and troughs in the amplitude of the oscillations as a function of 
wavenumber. Additionally the inflationary prediction of an approximately scale- 
invariant amplitude of the initial perturbations implies roughly scale-invariant 
oscillations in the power spectrum. And inflation generically predicts a flat Uni- 
verse. These are all falsifiablc i)rcdictions of the simplest inflationary models and 
they have withstood the test against observations to date. 

The energy contents of the Universe before recombination all leave their distinct 
signatures on the oscillations as discussed in §3.3-§3.5. In particular, the cold dark 
matter and baryon signatures have now been seen in the data (Halverson et al 
2001; Netterficld et al 2001; Lee et al 2001). The coupling l)ctwccu electrons and 
photons is not perfect, especially as one approaches the epoch of recombination. 
As discussed in §3.6, this imperfect coupling leads to damping in the anisotropy 
spectrum: very small scale inhomogeneities are smoothed out. The damping 
phenomenon has now been observed by the CBI experiment (Padin et al 2001). 
Importantly, fluid imperfections also generate linear polarization as covered in 
§3.7. Because the imperfection is minimal and appears only at small scales, the 
polarization generated is small and has not been detected to date. 
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After recombination the photons basically travel freely to us today, so the 
problem of translating the acoustic inhomogeneities in the photon distribution 
at recombination to the anisotropy spectrum today is simply one of projection. 
This projection depends almost completely on one number, the angular diameter 
distance between us and the surface of last scattering. That number depends on 
the energy contents of the Universe after recombination through the expansion 
rate. The hand waving projection argument of §3.1 is formalized in §3.8, in the 
process introducing the popular code used to compute anisotropics, cmbfast. 
Finally, we discuss the sensitivity of the acoustic peaks to cosmological parameters 
in §3.9. 

3. 1 Basics 

For pedagogical purposes, let us begin with an idealization of a perfect photon- 
baryon fluid and neglect the dynamical effects of gravity and the baryons. Per- 
turbations in this perfect fluid can be described by a simple continuity and an 
Euler equation that encapsulate the basic properties of acoustic oscillations. 

The discussion of acoustic oscillations will take place exclusively in Fourier 
space. For example, we decompose the monopole of the temperature fleld into 

^e^''--e(k), (6) 

and omit the subscript oo on the Fourier amplitude. Since perturbations are very 
small, the evolution equations are linear, and different Fourier modes evolve in- 
dependently. Therefore, instead of partial differential equations for a field 0(x), 
we have ordinary differential equations for 6(k). In fact, due to rotational sym- 
metry, all 0(k) for a given k obey the same equations. Here and in the following 
sections, we omit the wavenumber argument k where no confusion with physical 
space quantities will arise. 

Temperature perturbations in Fourier space obey 

e = -jK> (7) 

This equation for the photon temperature G, which does indeed look like the 
familiar continuity equation in Fourier space (derivatives V become wavenumbers 
fe), has a number of subtleties hidden in it, due to the cosmological setting. 

First, the "time" derivative here is actually with respect to conformal time -q = 
f dt/a{t). Since we are working in units in which the speed of light c = 1, i] is 
also the maximum comoving distance a particle could have traveled since t = 0. 
It is often called the comoving horizon or more specifically the comoving particle 
horizon. The physical horizon is a times the comoving horizon. 

Second, the photon fluid velocity here v^^ has been written as a scalar instead 
of a vector. In the early universe, only the velocity component parallel to the 
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wavevector k is expected to be important, since they alone have a source in 
gravity. Specifically, = — w-yk. In terms of the moments introduced in §2, v-y 
represents a dipole moment directed along k. The factor of 1/3 comes about since 
continuity conserves photon number not temperature and the number density 
rin/ oc T'\ Finally, we emphasize that, for the time being, we are neglecting the 
effects of gravity. 

The Euler equation for a fluid is an expression of momentum conservation. 
The momentum density of the photons is {pj + Pj)Vj, where the photon pressure 
Pj = Pj/3. In the absence of gravity and viscous fluid imperfections, pressure 
gradients Vp-y = Vp-y/3 supply the only force. Since Pj oc T^, this becomes 
4A:0p^/3 in Fourier space. The Euler equation then becomes 

v-y = ke . (8) 

Differentiating the continuity equation and inserting the Euler equation yields 
the most basic form of the oscillator equation 

e + c^k^e = 0, (9) 

where Cs = \/p/ p = l/\/3 is the sound speed in the (dynamically baryon-free) 
fluid. What this equation says is that pressure gradients act as a restoring force 
to any initial perturbation in the system which thereafter oscillate at the speed 
of sound. Physically these temperature oscillations represent the heating and 
cooling of a fluid that is compressed and rarefied by a standing sound or acoustic 
wave. This behavior continues until recombination. Assuming negligible initial 
velocity perturbations, we have a temperature distribution at recombination of 

e(r?,) = e(0)cos(fcs,), (10) 

where s = J Cgdr] ~ ^//VS is the distance sound can travel by r;, usually called 
the sound horizon. Asterisks denote evaluation at recombination a^. 

In the limit of scales large compared with the sound horizon /cs ^ 1, the 
perturbation is frozen into its initial conditions. This is the gist of the statement 
that the large-scale anisotropies measured by COBE directly measure the initial 
conditions. On small scales, the amplitude of the Fourier modes will exhibit 
temporal oscillations, as shown in Figure 1 [with 4" = 0, = 30(0) for this 
idealization]. Modes that are caught at maxima or minima of their oscillation 
at recombination correspond to peaks in the power, i.e. the variance of 0(fc, rj^,). 
Because sound takes half as long to travel half as far, modes corresponding to 
peaks follow a harmonic relationship fc„ = nn/s^, where n is an integer (see 
Figure la). 

How does this spectrum of inhomogeneities at recombination appear to us 
today? Roughly speaking, a spatial inhomogeneity in the CMB temperature of 
wavelength A appears as an angular anisotropy of scale 6 ~ X/D where D(z) is 
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s/st. s/s,. 

Figure 1: Idealized acoustic oscillations, (a) Peak scales: the wavemode that completes half 
an oscillation by recombination sets the physical scale of the first peak. Both minima and 
maxima correspond to peaks in power (dashed lines, absolute value) and so higher peaks are 
integral multiples of this scale with equal height. Plotted here is the idealization of Equation 
(15) (constant potentials, no baryon loading), (b) Baryon loading. Baryon loading boosts the 
amplitudes of every other oscillation. Plotted here is the idealization of Equation (16) (constant 
potentials and baryon loading i? — 1/6) for the third peak. 



the comoving angular diameter distance from the observer to redshift z. We will 
address this issue more formally in §3.8. In a fiat universe, Dt = ?7o ^ V* ~ Vo, 
where r]Q = r){z = 0). In harmonic space, the relationship implies a coherent 
series of acoustic peaks in the anisotropy spectrum, located at 

in ^ nl^, £a = nD^/s-, . (11) 

To get a feel for where these features should appear, note that in a flat matter 
dominated universe t] oc {1 + z)~^/^ so that "q^/riQ 1/30 ~ 2°. Equivalently 
l\ ~ 200. Notice that since we are measuring ratios of distances the absolute 
distance scale drops out; we shall see in §3.5 that the Hubble constant sneaks 
back into the problem because the Universe is not fully matter-dominated at 
recombination. 

In a spatially curved universe, the angular diameter distance no longer equals 
the coordinate distance making the peak locations sensitive to the spatial curva- 
ture of the Universe (Doroshkevich et al 1978; Kamionkowski et al 1994). Con- 
sider first a closed universe with radius of curvature R = HQ^\Q,tot — 1|~^/^. Sup- 
pressing one spatial coordinate yields a 2-sphere geometry with the observer situ- 
ated at the pole (see Figure 2). Light travels on lines of longitude. A physical scale 
A at fixed latitude given by the polar angle 6 subtends an angle a = X/RsinO. 
For a <C 1, a Euclidean analysis would infer a distance D = RsinO, even though 
the coordinate distance along the arc is d = dR; thus 



D = Rsm{d/R) 



(12) 
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Figure 2: Angular diameter distance. In a closed universe, objects are further than they appear 
to be from Euclidean (flat) expectations corresponding to the difference between coordinate 
distance d and angular diameter distance D. Consequently, at a fixed coordinate distance, a 
given angle corresponds to a smaller spatial scale in a closed universe. Acoustic peaks therefore 
appear at larger angles or lovi'er ^ in a closed universe. The converse is true for an open universe. 



For open universes, simply replace sin with sinh. The result is that objects in an 
open (closed) universe are closer (further) than they appear, as if seen through a 
lens. In fact one way of viewing this effect is as the gravitational lensing due to 
the background density (c.f. §4.2.4). A given comoving scale at a fixed distance 
subtends a larger (smaller) angle in a closed (open) universe than a flat universe. 
This strong scaling with spatial curvature indicates that the observed first peak 
at £i ~ 200 constrains the geometry to be nearly spatially flat. We will implicitly 
assume spatial flatness in the following sections unless otherwise stated. 

Finally in a flat dark energy dominated universe, the conformal age of the 
Universe decreases approximately as r;o ^ )7o(l + InH^''*^). For reasonable Qm, 
this causes only a small shift of £i to lower multipoles (see Plate 4) relative to the 
effect of curvature. Combined with the effect of the radiation near recombination, 
the peak locations provides a means to measure the physical age to of a flat 
universe (Hu et al 2001). 

3.2 Initial Conditions 

As suggested above, observations of the location of the first peak strongly point 
to a flat universe. This is encouraging news for adherents of inflation, a theory 
which initially predicted Htot = 1 at a time when few astronomers would sign 
on to such a high value (see Liddle & Lyth 1993 for a review). However, the 
argument for inflation goes beyond the confirmation of flatness. In particular, 
the discussion of the last subsection begs the question: whence B(0), the initial 
conditions of the temperature fluctuations? The answer requires the inclusion of 
gravity and considerations of causality which point to inflation as the origin of 
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structure in the Universe. 

The calculations of the typical angular scale of the acoustic oscillations in the 
last section are familiar in another context: the horizon problem. Because the 
sound speed is near the speed of light, the degree scale also marks the extent of 
a causally connected region or particle horizon at recombination. For the picture 
in the last section to hold, the perturbations must have been laid down while 
the scales in question were still far outside the particle horizon^. The recent 
observational verification of this basic peak structure presents a problem poten- 
tially more serious than the original horizon problem of approximate isotropy: 
the mechanism which smooths fluctuations in the Universe must also regenerate 
them with superhorizon sized correlations at the 10~^ level. Inflation is an idea 
that solves both problems simultaneously. 

The inflationary paradigm postulates that an early phase of near exponential 
expansion of the Universe was driven by a form of energy with negative pressure. 
In most models, this energy is usually provided by the potential energy of a scalar 
field. The inflationary era brings the observable universe to a nearly smooth and 
spatially flat state. Nonetheless, quantum fluctuations in the scalar field are un- 
avoidable and also carried to large physical scales by the expansion. Because an 
exponential expansion is self-similar in time, the fluctuations are scale-invariant, 
i.e. in each logarithmic interval in scale the contribution to the variance of the 
fluctuations is equal. Since the scalar field carries the energy density of the Uni- 
verse during inflation, its fluctuations induce variations in the spatial curvature 
(Guth & Pi 1985; Hawking 1982; Bardeen et al 1983). Instead of perfect flatness, 
inflation predicts that each scale will resemble a very slightly open or closed uni- 
verse. This fluctuation in the geometry of the Universe is essentially frozen in 
while the perturbation is outside the horizon (Bardeen 1980). 

Formally, curvature fluctuations are perturbations to the space-space piece 
of the metric. In a Newtonian coordinate system, or gauge, where the metric 
is diagonal, the spatial curvature fluctuation is called 6gij = 2a^<I>5y (see e.g. 
Ma & Bertschinger 1995). The more familiar Newtonian potential is the time- 
time fluctuation 6gtt = 2^' and is approximately ~ — <I>. Approximate scale 
invariance then says that A| = P^{k)/2Tv^ oc k"~^ where P$(fc) is the power 
spectrum of <1> and the tilt n ~ 1. 

Now let us relate the inflationary prediction of scale-invariant curvature fluctu- 
ations to the initial temperature fluctuations. Newtonian intuition based on the 
Poisson equation fc^<i> = AwGa^Sp tells us that on large scales (small k) density 
and hence temperature fluctuations should be negligible compared with Newto- 
nian potential. General relativity says otherwise because the Newtonian potential 
is also a time-time fluctuation in the metric. It corresponds to a temporal shift of 
5t/t = 9. The CMB temperature varies as the inverse of the scale factor, which 

^Recall that the comoving scale k does not vary with time. At very early times, then, the 
wavelengh is much larger than the horizon r). 
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in turn depends on time as a (x j2/[3(i+p/p)] Therefore, the fractional change in 
the CMB temperature 

e = -^ = -2fi + Er^^. (13) 

a 3\ pj t ^ ^ 

Thus, a temporal shift produces a temperature perturbation of — ^'/2 in the ra- 
diation dominated era (when p = p/3) and — 25'/3 in the matter dominated 
epoch (p = 0) (Peacock 1991; White & Hu 1997). The initial temperature per- 
turbation is therefore inextricably hnked with the initial gravitational potential 
perturbation. Inflation predicts scale-invariant initial fluctuations in both the 
CMB temperature and the spatial curvature in the Newtonian gauge. 

Alternate models which seek to obey the causality can generate curvature fluc- 
tuations only inside the particle horizon. Because the perturbations are then not 
generated at the same epoch independent of scale, there is no longer a unique 
relationship between the phase of the oscillators. That is, the argument of the 
cosine in Equation (10) becomes ks^, + 0(k), where is a phase which can in 
principle be different for different wavevectors, even those with the same mag- 
nitude k. This can lead to temporal incoherence in the oscillations and hence a 
washing out of the acoustic peaks (Albrecht et al 1996), most notably in cosmo- 
logical defect models (Allen et al 1997; Seljak et al 1997). Complete incoherence 
is not a strict requirement of causality since there arc other ways to synch up 
the oscillations. For example, many isocurvature models, where the initial spa- 
tial curvature is unperturbed, are coherent since their oscillations begin with the 
generation of curvature fluctuations at horizon crossing (Hu & White 1996). Still 
they typically have 4> ^ (c.f. Turok 1996). Independent of the angular diameter 
distance D*, the ratio of the peak locations gives the phase: f i : : ^3 1 : 2 : 3 
for (j) = 0. Likewise independent of a constant phase, the spacing of the peaks 
— (-n-i = (-A gives a measure of the angular diameter distance (Hu & White 
1996). The observations, which indicate coherent oscillations with <(> = 0, there- 
fore have provided a non-trivial test of the inflationary paradigm and supplied 
a substantially more stringent version of the horizon problem for contenders to 
solve. 



3.3 Gravitational Forcing 

We saw above that fluctuations in a scalar field during inflation get turned into 
temperature fluctuations via the intermediary of gravity. Gravity affects in 
more ways than this. The Newtonian potential and spatial curvature alter the 
acoustic oscillations by providing a gravitational force on the oscillator. The 
Euler equation (8) gains a term on the rhs due to the gradient of the potential 
k'^. The main effect of gravity then is to make the oscillations a competition 
between pressure gradients kQ and potential gradients k^ with an equilibrium 
when -I- = 0. 
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Gravity also changes the continuity equation. Since the Newtonian curvature 
is essentially a perturbation to the scale factor, changes in its value also generate 
temperature perturbations by analogy to the cosmological redshift 56 = —5$ 
and so the continuity equation (7) gains a contribution of — $ on the rhs. 

These two effects bring the oscillator equation (9) to 

1.2 

e-hc^fc^e = (14) 

In a flat universe and in the absence of pressure, $ and are constant. Also, 
in the absence of baryons, = 1/3 so the new oscillator equation is identical to 
Equation (9) with O replaced by 0-1-5'. The solution in the matter dominated 
epoch is then 

[e + *](?7) = [e + *](j],nd)cOs(fcs) 

= ^*(r7md) cos(fcs) • (15) 

where rj^d represents the start of the matter dominated epoch (see Figure la). We 
have used the matter dominated "initial conditions" for given in the previous 
section assuming large scales, feSmd ^ 1- 

The results from the idealization of §3.1 carry through with a few exceptions. 
Even without an initial temperature fluctuation to displace the oscillator, acoustic 
oscillations would arise by the infall and compression of the fluid into gravita- 
tional potential wells. Since it is the effective temperature -I- ^ that oscillates, 
they occur even if 0(0) = 0. The quantity -I- 5" can be thought of as an effec- 
tive temperature in another way: after recombination, photons must climb out 
of the potential well to the observer and thus suffer a gravitational redshift of 
AT/T = ^t. The effective temperature fluctuation is therefore also the observed 
temperature fluctuation. We now see that the large scale limit of Equation (15) 
recovers the famous Sachs- Wolfe result that the observed temperature perturba- 
tion is ^/3 and overdense regions correspond to cold spots on the sky (Sachs 
& Wolfe 1967). When 5' < 0, although is positive, the effective temperature 
-I- is negative. The plasma begins effectively rarefied in gravitational po- 
tential wells. As gravity compresses the fluid and pressure resists, rarefaction 
becomes compression and rarefaction again. The first peak corresponds to the 
mode that is caught in its first compression by recombination. The second peak 
at roughly half the wavelength corresponds to the mode that went through a full 
cycle of compression and rarefaction by recombination. We will use this language 
of the compression and rarefaction phase inside initially overdense regions but 
one should bear in mind that there are an equal number of initially underdense 
regions with the opposite phase. 
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3.4 Baryon Loading 

So far we have been neglecting the baryons in the dynamics of the acoustic oscil- 
lations. To see whether this is a reasonable approximation consider the photon- 
baryon momentum density ratio R = {pb + Pb) / {p-y + P-y) ~ 30Qbh^{z/W^)~^ . For 
typical values of the baryon density this number is of order unity at recombina- 
tion and so we expect baryonic effects to begin appearing in the oscillations just 
as they are frozen in. 

Baryons are conceptually easy to include in the evolution equations since their 
momentum density provides extra inertia in the joint Euler equation for pressure 
and potential gradients to overcome. Since inertial and gravitational mass are 
equal, all terms in the Euler equation save the pressure gradient are multiplied 
by 1 + -R leading to the revised oscillator equation (Hu & Sugiyama 1995) 

c?^(cr^e) + c?fc^e = -^*-c?J-(cj2<i>), (le) 

where we have used the fact that the sound speed is reduced by the baryons to 
Cs = l/V3(l + i?). 

To get a feel for the implications of the baryons take the limit of constant R, 
$ and Then cP{R9)/dri'^(= 0) may be added to the left hand side to again 
put the oscillator equation in the form of Equation (9) with G — > B + (1 + 
The solution then becomes 

[Q + {1 + R)-f]{ri) = [e + (l + i?)*]('/y^d)cos(M- (17) 

Aside from the lowering of the sound speed which decreases the sound horizon, 
baryons have two distinguishing effects: they enhance the amplitude of the oscil- 
lations and shift the equilibrium point to = —(1 + i?)*I' (see Figure lb). These 
two effects are intimately related and are easy to understand since the equations 
are exactly those of a mass m = 1 + i? on a spring in a constant gravitational 
field. For the same initial conditions, increasing the mass causes the oscillator to 
fall further in the gravitational field leading to larger oscillations and a shifted 
zero point. 

The shifting of the zero point of the oscillator has significant phenomenological 
consequences. Since it is still the effective temperature O + 'I' that is the observed 
temperature, the zero point shift breaks the symmetry of the oscillations. The 
baryons enhance only the compressional phase, i.e. every other peak. For the 
working cosmological model these are the first, third, fifth... Physically, the 
extra gravity provided by the baryons enhance compression into potential wells. 

These qualitative results remain true in the presence of a time-variable R. 
An additional effect arises due to the adiabatic damping of an oscillator with a 
time- variable mass. Since the energy/frequency of an oscillator is an adiabatic 
invariant, the amplitude must decay as (1 + 7?)"^/*. This can also be understood 
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Figure 3: Radiation driving and diffusion damping. The decay of the potential ^ drives the 
oscillator in the radiation dominated epoch. Diffusion generates viscosity ttj, i.e. a quadrupole 
moment in the temperature, which damps oscillations and generates polarization. Plotted here 
is the numerical solution to Equation (18) and Equation (19) for a mode with wavelength much 
smaller than the sound horizon at decoupling, ks* ^ 1. 

by expanding the time derivatives in Equation (16) and identifying the RQ term 
as the remnant of the familiar expansion drag on baryon velocities. 

3.5 Radiation Driving 

We have hitherto also been neglecting the energy density of the radiation in 
comparison to the matter. The matter-to-radiation ratio scales as Pm/pr ~ 
24Q„/i2(2/l03)-i and so is also of order unity at recombination for reasonable 
parameters. Moreover fiuctuations corresponding to the higher peaks entered the 
sound horizon at an earlier time, during radiation domination. 

Including the radiation changes the expansion rate of the Universe and hence 
the physical scale of the sound horizon at recombination. It introduces yet an- 
other potential ambiguity in the interpretation of the location of the peaks. For- 
tunately, the matter-radiation ratio has another effect in the power spectrum by 
which it can be distinguished. Radiation drives the acoustic oscillations by mak- 
ing the gravitational force evolve with time (Hu & Sugiyama 1995). Matter does 
not. 

The exact evolution of the potentials is determined by the relativistic Poisson 
equation. But qualitatively, we know that the background density is decreasing 
with time, so unless the density fluctuations in the dominant component grow 
unimpeded by pressure, potentials will decay. In particular, in the radiation 
dominated era once pressure begins to fight gravity at the first compressional 
maxima of the wave, the Newtonian gravitational potential and spatial curvature 
must decay (see Figure 3). 
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This decay actually drives the oscillations; it is timed to leave the fluid max- 
imally compressed with no gravitational potential to fight as it turns around. 
The net efl^ect is doubled since the redshifting from the spatial metric fluctua- 
tion $ also goes away at the same time. When the Universe becomes matter 
dominated the gravitational potential is no longer determined by photon-baryon 
density perturbations but by the prcssureless cold dark matter. Therefore, the 
amplitudes of the acoustic peaks increase as the cold dark matter-to-radiation 
ratio decreases (Seljak 1994; Hu & Sugiyama 1995). Density perturbations in 
any form of radiation will stop growing around horizon crossing and lead to this 
effect. The net result is that across the horizon scale at matter radiation equality 
(fceq = (4 — 2\/2)/r}eq) tlic acoustic amplitude increases by a factor of 4-5 (Hu & 
Sugiyama 1996). By eliminating gravitational potentials, photon-baryon acous- 
tic oscillations eliminate the alternating peak heights from baryon loading. The 
observed high third peak (Halverson et al 2001) is a good indication that cold 
dark matter both exists and dominates the energy density at recombination. 

3. 6 Damping 

The photon-baryon fluid has slight imperfections corresponding to shear viscosity 
and heat conduction in the fluid (Weinberg 1971). These imperfections damp 
acoustic oscillations. To consider these effects, we now present the equations of 
motion of the system in their full form, including separate continuity and Euler 
equations for the baryons. Formally the continuity and Euler equations follow 
from the covariant conservation of the joint stress-energy tensor of the photon- 
baryon fluid. Because photon and baryon numbers are separately conserved, the 
continuity equations are unchanged, 

e = ~v^-^, 5b = -kvb-3^, (18) 

where Sb and Vb are the density perturbation and fluid velocity of the baryons. 
The Euler equations contain qualitatively new terms 

k 

Vb = ~vb + k^ + f{vy-Vb)/R. (19) 

For the baryons the first term on the right accounts for cosmological expansion, 
which makes momenta decay as a~^. The third term on the right accounts for 
momentum exchange in the Thomson scattering between photons and electrons 
(protons are very tightly coupled to electrons via Coulomb scattering), with f = 
UeO-TCi the differential Thomson optical depth, and is compensated by its opposite 
in the photon Euler equation. These terms are the origin of heat conduction 
imperfections. If the medium is optically thick across a wavelength, f/k ^ 1 and 



20 



Hu & DodelsoD 



the photons and baryons cannot slip past each other. As it becomes optically 
thin, slippage dissipates the fluctuations. 

In the photon Euler equation there is an extra force on the rhs due to anisotropic 
stress gradients or radiation viscosity in the fluid, n^. The anisotropic stress is 
directly proportional to the quadrupole moment of the photon temperature dis- 
tribution. A quadrupole moment is established by gradients in v^/ as photons 
from say neighboring temperature crests meet at a trough (see Plate 3, inset). 
However it is destroyed by scattering. Thus ir-y = 2(kv-y/f)Ay, where the order 
unity constant can be derived from the Boltzmann equation = 16/15 (Kaiser 
1983). Its evolution is shown in Figure 3. With the continuity Equation (7), 
kvj « —30 and so viscosity takes the form of a damping term. The heat conduc- 
tion term can be shown to have a similar eft'cct by expanding the Euler equations 
in k/f. The flnal oscillator equation including both terms becomes 

c',^{c7'e) + ^[A^ + A^]@ + cy@ = -^vP-c^A(cr2$), (20) 

where the heat conduction coefficient Ah = R^/{1 + R). Thus we expect the 
inhomogeneities to be damped by a exponential factor of order e~^^'^/'^ (see Figure 
3). The damping scale fed is thus of order ^f/rj, corresponding to the geometric 
mean of the horizon and the mean free path. Damping can be thought of as the 
result of the random walk in the baryons that takes photons from hot regions into 
cold and vice- versa (Silk 1968). Detailed numerical integration of the equations of 
motion are required to track the rapid growth of the mean free path and damping 
length through recombination itself. These calculations show that the damping 
scale is of order fcas* « 10 leading to a substantial suppression of the oscillations 
beyond the third peak. 

How does this suppression depend on the cosmological parameters? As the 
matter density Qmh'^ increases, the horizon 77* decreases since the expansion rate 
goes up. Since the diffusion length is proportional to it too decreases as 
the matter density goes up but not as much as the angular diameter distance 
which is also inversely proportional to the expansion rate. Thus, more matter 
translates into more damping at a fixed multipole moment; conversely, it cor- 
responds to slightly less damping at a fixed peak number. The dependence on 
baryons is controlled by the mean free path which is in turn controlled by the free 
electron density: the increase in electron density due to an increase in the baryons 
is partially offset by a decrease in the ionization fraction due to recombination. 
The net result under the Saha approximation is that the damping length scales 
approximately as (ft^/i^)"^/^. Accurate fitting formulae for this scale in terms of 
cosmological parameters can be found in (Hu & White 1997c). 
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3. 7 Polarization 

The dissipation of the acoustic osciUations leaves a signature in the polarization 
of CMB in its wake (see e.g. Hu & White 1997a and references therein for a more 
complete treatment). Much like reflection off of a surface, Thomson scattering 
induces a linear polarization in the scattered radiation. Consider incoming ra- 
diation in the — x direction scattered at right angles into the z direction (see 
Plate 2, left panel). Heuristically, incoming radiation shakes an electron in the 
direction of its electric field vector or polarization e' causing it to radiate with 
an outgoing polarization parallel to that direction. However since the outgoing 
polarization e must be orthogonal to the outgoing direction, incoming radiation 
that is polarized parallel to the outgoing direction cannot scatter leaving only 
one polarization state. More generally, the Thomson differential cross section 

Unlike the reflection of sunlight off of a surface, the incoming radiation comes 
from all angles. If it were completely isotropic in intensity, radiation coming along 

the y would provide the polarization state that is missing from that coming along 
X leaving the net outgoing radiation unpolarized. Only a quadrupole temperature 
anisotropy in the radiation generates a net linear polarization from Thomson 
scattering. As we have seen, a quadrupole can only be generated causally by the 
motion of photons and then only if the Universe is optically thin to Thomson 
scattering across this scale (i.e. it is inversely proportional to f). Polarization 
generation suffers from a Catch-22: the scattering which generates polarization 
also suppresses its quadrupole source. 

The fact that the polarization strength is of order the quadrupole explains the 
shape and height of the polarization spectra in Plate lb. The monopole and 
dipole and are of the same order of magnitude at recombination, but their 
oscillations are n/2 out of phase as follows from Equation (9) and Equation (10). 
Since the quadrupole is of order kv^/f (see Figure 3), the polarization spectrum 
should be smaller than the temperature spectrum by a factor of order fc/f at 
recombination. As in the case of the damping, the precise value requires numerical 
work (Bond & Efstatliiou 1987) since f changes so rapidly near recombination. 
Calculations show a steady rise in the polarized fraction with increasing Z or fe to a 
maximum of about ten percent before damping destroys the oscillations and hence 
the dipole source. Since v-y is out of phase with the monopole, the polarization 
peaks should also be out of phase with the temperature peaks. Indeed, Plate lb 
shows that this is the case. Furthermore, the phase relation also tells us that the 
polarization is correlated with the temperature perturbations. The correlation 
power C®^ being the product of the two, exhibits oscillations at twice the acoustic 
frequency. 

Until now, we have focused on the polarization strength without regard to its 

orientation. The orientation, like a 2 dimensional vector, is described by two 
components E and B. The E and B decomposition is simplest to visualize in 
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the small scale limit, where spherical harmonic analysis coincides with Fourier 
analysis (Seljak 1997). Then the wavevector k picks out a preferred direction 
against which the polarization direction is measured (see Plate 2, right panel). 
Since the linear polarization is a "headless vector" that remains unchanged upon 
a 180° rotation, the two numbers E and B that define it represent polarization 
aligned or orthogonal with the wavevector (positive and negative E) and crossed 
at ±45° (positive and negative B). 

In linear theory, scalar perturbations like the gravitational potential and tem- 
perature perturbations have only one intrinsic direction associated with them, 
that provided by k, and the orientation of the polarization inevitably takes it 
cue from that one direction, thereby producing an i?— mode. The generalization 
to an all-sky characterization of the polarization changes none of these qualita- 
tive features. The E—mode and the B— mode are formally distinguished by the 
orientation of the Hessian of the Stokes parameters which define the direction of 
the polarization itself. This geometric distinction is preserved under summation 
of all Fourier modes as well as the generalization of Fourier analysis to spherical 
harmonic analysis. 

The acoustic peaks in the polarization appear exclusively in the EE power 
spectrum of Equation (5). This distinction is very useful as it allows a clean 

separation of this effect from those occuring beyond the scope of the linear per- 
turbation theory of scalar fluctuations: in particular, gravitational waves (see 
§4.2.3) and gravitational lensing (see §4.2.4). Moreover, in the working cosmo- 
logical model, the polarization peaks and correlation are precise predictions of the 
temperature peaks as they depend on the same physics. As such their detection 
would represent a sharp test on the implicit assumptions of the working model, 
especially its initial conditions and ionization history. 

3.8 Integral Approach 

The discussion in the previous sections suffices for a qualitative understanding 
of the acoustic peaks in the power spectra of the temperature and polarization 
anisotropics. To refine this treatment we must consider more carefully the sources 
of anisotropics and their projection into multipole moments. 

Because the description of the acoustic oscillations takes place in Fourier space, 
the projection of inhomogeneities at recombination onto anisotropies today has 
an added level of complexity. An observer today sees the acoustic oscillations 
in offoctivo temperature as they appeared on a spherical shell at x = D^h at 
recombination, where n is the direction vector, and = rjQ — rj^, is the distance 
light can travel between recombination and the present (see Plate 3). Having 
solved for the Fourier amplitude [G-l-3'](fc, 77*), we can expand the exponential in 
Equation (6) in terms of spherical harmonics, so the observed anisotropy today 
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Plate 2: Polarization generation and classification. Left: Tlionison scattering of quadrupole 
temperature anisotropies (depicted here in the x — y plane) generates linear polarization. Right: 
Polarization in the x — y plane along the outgoing z axis. The component of the polarization 
that is parallel or perpendicular to the wavevector k is called the £^-mode and the one at 45° 
angles is called the B-mode. 
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where the projected source ae{k) = [0 + *I'](k, rj.f)ji(kDi,). Because the spherical 
harmonics are orthogonal, Equation (1) imphes that Qirn today is given by the 
integral in square brackets today. A given plane wave actually produces a range of 
anisotropies in angular scale as is obvious from Plate 3. The one-to-one mapping 
between wavenumber and multipole moment described in §3.1 is only approxi- 
mately true and comes from the fact that the spherical Bessel function jiikD^) is 
strongly peaked at fc_D* I. Notice that this peak corresponds to contributions 
in the direction orthogonal to the wavevector where the correspondence between 
I and k is one-to-one (see Plate 3). 

Projection is less straightforward for other sources of anisotropy. We have 
hitherto neglected the fact that the acoustic motion of the photon-baryon fluid 
also produces a Doppler shift in the radiation that appears to the observer as 
a temperature anisotropy as well. In fact, we argued above that ~ is 
of comparable magnitude but out of phase with the effective temperature. If 
the Doppler effect projected in the same way as the effective temperature, it 
would wash out the acoustic peaks. However, the Doppler effect has a directional 
dependence as well since it is only the line-of-sight velocity that produces the 
effect. Formally, it is a dipole source of temperature anisotropies and hence 
has an £ = 1 structure. The coupling of the dipole and plane wave angular 
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Plate 3: Integral approach. CMB anisotropies can be thought of as the line-of-sight projection 
of various sources of plane wave temperature and polarization fluctuations: the acoustic effective 
temperature and velocity or Doppler effect (see §3.8), the quadrupole sources of polarization (see 
§3.7) and secondary sources (see §4.2, §4.3). Secondary contributions differ in that the region 
over vi'hich they contribute is thick compared with the last scattering surface at recombination 
and the typical wavelength of a perturbation. 
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momenta imply that in the projection of the Doppler effect involves a combination 
of je±i that may be rewritten as j({x) = dje{x)/dx. The structure of j'^ lacks a 
strong peak at x = i. Physically this corresponds to the fact that the velocity 
is irrotational and hence has no component in the direction orthogonal to the 
wavevector (sec Plate 3). Correspondingly, the Doppler effect cannot produce 
strong peak structures (Hu & Sugiyama 1995). The observed peaks must be 
acoustic peaks in the effective temperature not "Doppler peaks" . 

There is one more subtlety involved when passing from acoustic oscillations to 
anisotropies. Recall from §3.5 that radiation leads to decay of the gravitational 
potentials. Residual radiation after decoupling therefore implies that the effective 
temperature is not precisely [0 + <I'](r;,). The photons actually have slightly 
shallower potentials to climb out of and lose the pcrturbative analogue of the 
cosmological redshift, so the [© + $](?7») overestimates the difference between the 
true photon temperature and the observed temperature. This effect of course is 
already in the continuity c<iuation for the monopole Equation (18) and so the 
source in Equation (21) gets generalized to 

ae{k) = [6 + *] {v*)ji(.kD,) + vb{k,-n,)j[{kD^) + / dr]{^ - ^)ji{kD). (22) 

■/»). 

The last term vanishes for constant gravitational potentials, but is non-zero if 
residual radiation driving exists, as it will in low Q-mh"^ models. Note that residual 
radiation driving is particularly important because it adds in phase with the 
monopole: the potentials vary in time only near recombination, so the Bessel 
function can be set to jiikD^,) and removed from the r] integral. This complication 
has the effect of decreasing the multipole value of the first peak t\ as the matter- 
radiation ratio at recombination decreases (Hu & Sugiyama 1995). Finally, we 
mention that time varying potentials can also play a role at very late times due to 
non-linearities or the importance of a cosmological constant for example. Those 
contributions, to be discussed more in §4.2.1, are sometimes referred to as late 
Integrated Sachs- Wolfe effects, and do not add coherently with [0 -I- ^](?7*). 

Putting these expressions together and squaring, we obtain the power spectrum 
of the acoustic oscillations 

Ce = - f ^k'4{k). (23) 

IT J K 

This formulation of the anisotropies in terms of projections of sources with 
specific local angular structure can be completed to include all types of sources 
of temperature and polarization anisotropies at any given epoch in time linear 
or non-linear: the monopole, dipole and quadrupole sources arising from density 
perturbations, vorticity and gravitational waves (Hu & White 1997b). In a curved 
geometry one replaces the spherical Bessel functions with ultraspherical Bessel 
functions (Abbott & Schaefer 1986; Hu et al 1998). Precision in the predictions 
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of the observables is then limited only by the precision in the prediction of the 
sources. This formulation is ideal for cases where the sources are governed by 
non-linear physics even though the CMB responds linearly as we shall see in §4. 

Perhaps more importantly, the widely-used cmbfast code (Seljak & Zaldar- 
riaga 1996) exploits these properties to calculate the anisotropies in linear per- 
turbation efficiently. It numerically solves for the smoothly-varying sources on 
a sparse grid in wavenumber, interpolating in the integrals for a handful of ^'s 
in the smoothly varying Q. It has largely replaced the original ground breaking 
codes (Wilson & Silk 1981; Bond & Efstathiou 1984; Vittorio & Silk 1984) based 
on tracking the rapid temporal oscillations of the multipole moments that simply 
reflect structure in the spherical Bessel functions themselves. 



3.9 Parameter Sensitivity 

The phenomenology of the acoustic peaks in the temperature and polarization is 
essentially described by 4 observables and the initial conditions (Hu et al 1997). 
These are the angular extents of the sound horizon £a = TrDt/st, the particle hori- 
zon at matter radiation equality i^q = k^qD^ and the damping scale = k^Ds^ 
as well as the value of the baryon-photon momentum density ratio £a sets the 
spacing between of the peaks; ^eq and £4 compete to determine their amplitude 
through radiation driving and diffusion damping. R^, sets the baryon loading 
and, along with the potential well depths set by £eq, fixes the modulation of the 
even and odd peak heights. The initial conditions set the phase, or equivalently 
the location of the first peak in units of tg,, and an overall tilt n in the power 
spectrum. 

In the model of Plate 1, these numbers are = 301 {Ci = 0.73£a), ieq = 149, 
= 1332, 7?* = 0.57 and n = 1 and in this family of models the parameter 
sensitivity is approximately (Hu et al 2001) 

^ « -0.24^+0.07^-0.17^-1.1^, 

A4q ^ ^^AQ^h'^ ^_^^AQf^ _^ArJtot ^^^-^ 

^eq ^jjih^ ^tot 

and AR^/Rt ki 1.0AQi,h'^/Q,iih'^. Current observations indicate that = 304±4, 
4q = 168 ± 15, 4 = 1392 ± 18, fl* = 0.60 ± 0.06, and n = 0.96 ± 0.04 (Knox 
et al 2001; see also Wang et al 2001; Pryke et al 2001; de Bernardis et al 2001), if 
gravitational waves contributions are subdominant and the reionization redshift 
is low as assumed in the working cosmological model (see §2.1). 

The acoustic peaks therefore contain three rulers for the angular diameter 
distance test for curvature, i.e. deviations from fitot = 1- However contrary to 
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100 - (a) Curvature 



(b) Dark Energy 




Plate 4: Sensitivity of tlie acoustic temperature spectrum to four fundamental cosmological 
parameters (a) the curvature as quantified by Qtot (b) the dark energy as quantified by the 
cosmological constant Oa {wa = —1) (c) the physical baryon density Q^h^ (d) the physical 
matter density Qmh'^; all varied around a fiducial model of Qtot — 1, — 0.65, Qbh^ — 0.02, 
n^h'^ = 0.147, n = 1, Zri = 0, = 0. 
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popular belief, any one of these alone is not a standard ruler whose absolute 
scale is known even in the working cosmological model. This is refiected in the 
sensitivity of these scales to other cosmological parameters. For example, the 
dependence of £b. on and hence the Hubble constant is quite strong. But 

in combination with a measurement of the matter-radiation ratio from i^q, this 
degeneracy is broken. 

The weaker degeneracy of on the baryons can likewise be broken from a 
measurement of the baryon-plioton ratio . The damping scale provides an 
additional consistency check on the implicit assumptions in the working model, 
e.g. recombination and the energy contents of the Universe during this epoch. 
What makes the peaks so valuable for this test is that the rulers are standardize- 
able and contain a built-in consistency check. 

There remains a weak but perfect degeneracy between Htot and SIa because 
they both appear only in D*. This is called the angular diameter distance degen- 
eracy in the literature and can readily be generalized to dark energy components 
beyond the cosmological constant assumed here. Since the effect of Qa is in- 
trinsically so small, it only creates a correspondingly small ambiguity in f2tot for 
reasonable values of Q\. The down side is that dark energy can never be isolated 
through the peaks alone since it only takes a small amount of curvature to mimic 
its effects. The evidence for dark energy through the CMB comes about by al- 
lowing for external information. The most important is the nearly overwhelming 
direct evidence for Q„i < 1 from local structures in the Universe. The second is 
the measurements of a relatively high Hubble constant h ~ 0.7; combined with a 
relatively low Qjnh^ that is preferred in the CMB data, it implies Qm < 1 but at 
low significance currently. 

The upshot is that precise measurements of the acoustic peaks yield precise de- 
terminations of four fundamental parameters of the working cosmological model: 
Q,mh? , D*, and n. More generally, the first three can be replaced by ig., ^eq, 
£d and i?* to extend these results to models where the underlying assumptions 
of the working model are violated. 

4 BEYOND THE PEAKS 

Once the acoustic peaks in the temperature and polarization power spectra have 
been scaled, the days of splendid isolation of cosmic microwave background the- 
ory, analysis and experiment will have ended. Beyond and beneath the peaks 
lies a wealth of information about the evolution of structure in the Universe and 
its origin in the early universe. As CMB photons traverse the large scale struc- 
ture of the Universe on their journey from the recombination epoch, they pick 
up secondary temperature and polarization anisotropies. These depend on the 
intervening dark matter, dark energy, baryonic gas density and temperature dis- 
tributions, and even the existence of primordial gravity waves, so the potential 
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payoff of their detection is enormous. The price for this extended reach is the 
loss of the abihty both to make precise predictions, due to uncertain and/or non- 
linear physics, and to make precise measurements, due to the cosmic variance of 
the primary anisotropies and the relatively greater importance of galactic and 
extragalactic foregrounds. 

We begin in §4.1 with a discussion of the matter power spectrum to set the 
framework for the discussion of secondary anisotropies. Secondaries can be di- 
vided into two classes: those due to gravitational effects and those induced by 
scattering off of electrons. The former are treated in §4.2 and the latter in §4.3. 
Secondary anisotropies are often non-Gaussian, so they show up not only in the 
power spectra of §2, but in higher point functions as well. We briefly discuss non- 
Gaussian statistics in §4.4. All of these topics are subjects of current research to 
which this review can only serve as introduction. 

4-1 Matter Power Spectrum 

The same Vialance between pressure and gravity that is respoiisiblo for acoustic 
oscillations determines the power spectrum of fluctuations in the non-relativistic 
matter. This relationship is often obscured by focussing on the density fluctua- 
tions in the pressureless cold dark matter itself and we so we will instead consider 
the matter power spectrum from the perspective of the Newtonian potential. 

4.1.1 Physical Description After recombination, without the pressure of 
the photons, the baryons simply fall into the Newtonian potential wells with the 
cold dark matter, an event usually referred to as the end of the Compton drag 
epoch. We claimed in §3.5 that above the horizon at matter-radiation equality the 
potentials are nearly constant. This follows from the dynamics: where pressure 
gradients are negligible, infall into some initial potential causes a potential flow of 
vtot ^ {kri)'^i [see Equation (19)] and causes density enhancements by continuity 
of 5tot ~ —{krOv^ot ~ ^{kri)'^^i. The Poisson equation says that the potential 
at this later time ^ ~ — (fc'?)~^'5tot ^ so that this rate of growth is exactly 
right to keep the potential constant. Formally, this Newtonian argument only 
applies in general relativity for a particular choice of coordinates (Bardeen 1980), 
but the rule of thumb is that if what is driving the expansion (including spatial 
curvature) can also cluster unimpeded by pressure, the gravitational potential 
will remain constant. 

Because the potential is constant in the matter dominated epoch, the large- 
scale observations of COBE set the overall amislitude of the potential power 
spectrum today. Translated into density, this is the well-known COBE normal- 
ization. It is usually expressed in terms of 5h, the matter density perturbation 
at the Hubble scale today. Since the observed temperature fluctuation is approx- 
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imately 3'/3, 

#^A|.i4, (25) 

where the second equality follows from the Poisson equation in a fully matter 
dominated universe with £1^ = 1. The observed GOBE fluctuation of Ay « 28^K 
(Smoot et al 1992) implies 5h ~ "2 x 10~^. For corrections for 0,^, < 1 where the 
potential decays because the dominant driver of the expansion cannot cluster, 
see Bunn & White (1997). 

On scales below the horizon at matter-radiation equality, we have seen in §3.5 
that pressure gradients from the acoustic oscillations themselves impede the clus- 
tering of the dominant component, i.e. the photons, and lead to decay in the 
potential. Dark matter density perturbations remain but grow only logarithmi- 
cally from their value at horizon crossing, which (just as for large scales) is ap- 
proximately the initial potential, 5m ~ —"^i- The potential for modes that have 
entered the horizon already will therefore be suppressed by $ oc —dm/k'^ ~ ^i/k"^ 
at matter domination (neglecting the logarithmic growth) again according to the 
Poisson equation. The ratio of ^ at late times to its initial value is called the 
transfer function. On large scales, then, the transfer function is close to one, 
while it falls o£F as on small scales. If the baryons fraction Pb/ pm is sub- 
stantial, baryons alter the transfer function in two ways. First their inability to 
cluster below the sound horizon causes further decay in the potential between 
matter-radiation equality and the end of the Compton drag epoch. Secondly the 
acoustic oscillations in the baryonic velocity field kinematically cause acoustic 
wiggles in the transfer function (Hu & Sugiyama 1996). These wiggles in the 
matter power spectrum are related to the acoustic peaks in the CMB spectrum 
like twins separated at birth and are actively being pursued by the largest galax;y 
surveys (Percival et al 2001). For fitting formulae for the transfer function that 
include these effects see Eisenstein & Hu (1998). 



4.1.2 COSMOLOGICAL Implications The combination of the GOBE normal- 
ization, the matter transfer function and the near scale-invariant initial spectrum 
of fluctuations tells us that by the present fluctuations in the cold dark matter 
or baryon density fields will have gone non-linear for all scales k > 10~^/iMpc~"'. 
It is a great triumph of the standard cosmological paradigm that there is just 
enough growth between ~ 10^ and « = to explain structures in the Universe 
across a wide range of scales. 

In particular, since this non-linear scale also corresponds to galaxy clusters 
and meastu'ements of their abundance yields a robust measure of the power near 
this scale for a given matter density Q^. The agreement between the GOBE 
normalization and the cluster abundance at low fi„i ^ 0.3 — 0.4 and the observed 
Hubble constant h = 0.72 ± 0.08 (Freedman et al 2001) was pointed out imme- 
diately following the COBE result (e.g. White et al 1993; Bartlett & Silk 1993) 
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and is one of the strongest pieces of evidence for the parameters in the working 
cosmological model (Ostriker & Steinhardt 1995; Krauss & Turner 1995). 

More generally, the comparison between large-scale structure and the CMB is 
important in that it breaks degeneracies between effects due to deviations from 
power law initial conditions and the dynamics of the matter and energy contents 
of the Universe. Any dynamical effect that reduces the amplitude of the matter 
power spectrum corresponds to a decay in the Newtonian potential that boosts 
the level of anisotropy (see §3.5 and §4.2.1). Massive neutrinos are a good example 
of physics that drives the matter power spectrum down and the CMB spectrum 
up. 

The combination is even more fruitful in the relationship between the acoustic 
peaks and the baryon wiggles in the matter power spectrum. Our knowledge of 
the physical distance between adjacent wiggles provides the ultimate standard 
candle for cosmology (Eisenstein et al 1998). For example, at very low z, the 
radial distance out to a galaxy is cz/Hq. The unit of distance is therefore 
Mpc, and a knowledge of the true physical distance corresponds to a determi- 
nation of h. At higher redshifts, the radial distance depends sensitively on the 
background cosmology (especially the dark energy), so a future measurement of 
baryonic wiggles at z ~ 1 say would be a powerful test of dark energy models. 
To a lesser extent, the shape of the transfer function, which mainly depends on 
the matter-radiation scale in h Mpc~^, i.e. ilmh, is another standard ruler (see 
e.g. Tegmark et al 2001 for a recent assessment), more heralded than the wiggles, 
but less robust due to degeneracy with other cosmological parameters. 

For scales corresponding to A; ^ 10~^h Mpc~^, density fluctuations are non- 
linear by the present. Numerical A^-body simulations show that the dark matter 
is bound up in a hierarchy of virialized structures or lialos (see Bertscliinger 1998 
for a review). The statistical properties of the dark matter and the dark matter 
halos have been extensively studied in the working cosmological model. Less 
certain are the properties of the baryonic gas. We shall see that both enter into 
the consideration of secondary CMB anisotropies. 

4-2 Gravitational Secondaries 

Gravitational secondaries arise from two sources: the differential redshift from 

time-variable metric perturbations (Sachs & Wolfe 1967) and gravitational lens- 
ing. There are many examples of the former, one of which we have already 
encountered in §3.8 in the context of potential decay in the radiation dominated 
era. Such gravitational potential effects are usually called the integrated Sachs- 
Wolfe (ISW) effect in linear perturbation theory (§4.2.1), the Rees-Sciama (§4.2.2) 
effect in the non-linear regime, and the gravitational wave effect for tensor per- 
turbations (§4.2.3). Gravitational waves and lensing also produce B-modes in 
the polarization (see §3.7) by which they may be distinguished from acoustic 
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polarization. 

4.2.1 ISW Effect As we have seen in the previous section, the potential 
on a given scale decays whenever the expansion is dominated by a component 
whose effective density is smooth on that scale. This occiirs at late times in an 
D,rn < 1 model at the end of matter domination and the onset dark energy (or 
spatial curvature) domination. If the potential decays between the time a photon 
falls into a potential well and when it climbs out it gets a boost in temperature 
of due to the differential gravitational redshift and —(5$ « d'^ due to an 
accompanying contraction of the wavelength (see §3.3). 

Potential decay due to residual radiation was introduced in §3.8, but that due 
to dark energy or curvature at late times induces much different changes in the 
anisotropy spectrum. What makes the dark energy or curvature contributions 
different from those due to radiation is the longer length of time over which the 
potentials decay, on order the Hubble time today. Residual radiation produces its 
effect quickly, so the distance over which photons feel the effect is much smaller 
than the wavelength of the potential fluctuation. Recall that this meant that 
ji{kD) in the integral in Equation (23) could be set to ji{kDt) and removed from 
the integral. The final effect then is proportional to ji{kD^) and adds in phase 
with the monopole. 

The ISW projection, indeed the projection of all secondaries, is much different 
(see Plate 3). Since the duration of the potential change is much longer, pho- 
tons typically travel through many peaks and troughs of the perturbation. This 
cancellation implies that many modes have virtually no impact on the photon 
temperature. The only modes which do have an impact are those with wavevec- 
tors perpendicular to the line of sight, so that along the line of sight the photon 
does not pass through crests and troughs. What fraction of the modes contribute 
to the effect then? For a given wavenumber k and line of sight instead of the 
full spherical shell at radius Airk'^dk, only the ring 2nkdk with k _L n participate. 
Thus, the anisotropy induced is suppressed by a factor of k (or I = kD in angular 
space). Mathematically, this arises in the line-of-sight integral of Equation (23) 
from the integral over the oscillatory Bessel function / dxje{x) « (7r/2^)^/^ (see 
also Plate 3). 

The ISW effect thus generically shows up only at the lowest ^'s in the power 
spectrum (Kofman & Starobinskii 1985). This spectrum is shown in Plate 5. 
Secondary anisotropy predictions in this figure are for a model with fitot = li 
= 2/3, Qfjh'^ = 0.02, Q.mh'^ = 0.16, n = 1, Zn = 7 and inflationary energy 
scale Ei <^ IQ^^ GeV. The ISW effect is especially important in that it is ex- 
tremely sensitive to the dark energy: its amount, equation of state and clustering 
properties (Coble et al 1997; Caldwell et al 1998; Hu 1998). Unfortunately, being 
confined to the low multipoles, the ISW effect suffers severely from the cosmic 
variance in Equation (4) in its detectability. Perhaps more promising is its cor- 
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Plate 5: Secondary anisotropies. (a) Gravitational secondaries: ISW, lensing and Rees-Sciama 
(moving halo) effects, (b) Scattering secondaries: Doppler, density (S) and ionization (i) mod- 
ulated Doppler, and the SZ effects. Curves and model are described in the text. 

This type of cancellation behavior and corresponding suppression of small scale 
fluctuations is a common feature of secondary temperature and polarization an- 
isotropies from large-scale structure and is quantified by the Limber equation 
(Limber 1954) and its CMB generalization (Hu & White 1996; Hu 2000a). It 
is the central reason why secondary anisotropies tend to be smaller than the 
primary ones from 2,^ 10^ despite the intervening growth of structure. 

4.2.2 Rees-Sciama and Moving Halo Effects The ISW effect is linear 
in the perturbations. Cancellation of the ISW effect on small scales leaves second 
order and non-linear analogues in its wake (Rees & Sciama 1968). From a single 
isolated structure, the potential along the line of sight can change not only from 
evolution in the density profile but more importantly from its bulk motion across 
the line of sight. In the context of clusters of galaxies, this is called the moving 
cluster effect (Birkinshaw & Gull 1983). More generally, the bulk motion of dark 
matter halos of all masses contribute to this effect (Tuluie & Laguna 1995; Seljak 
1996b), and their clustering gives rise to a low level of anisotropies on a range of 
scales but is never the leading source of secondary anisotropies on any scale (see 
Plate 5a). 

4.2.3 Gravitational Waves A time- variable tensor metric perturbation 
similarly leaves an imprint in the temperature anisotropy (Sachs & Wolfe 1967). 
A tensor metric perturbation can be viewed as a standing gravitational wave and 
produces a quadrupolar distortion in the spatial metric. If its amplitude changes, 
it leaves a quadrupolar distortion in the CMB temperatui'e distribution (Polnarev 
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Figure 4: Gravitational waves and the energy scale of inflation Ei. Left: temperature and 
polarization spectra from an initial scale invariant gravitational wave spectrum with power 
oc Ef — (4 X lO^^GeV)'^. Right: 95% confidence upper limits statistically achievable on Ei and 
the scalar tilt n by the MAP and Planck satellites as well as an ideal experiment out to ^ — 3000 
in the presence of gravitational lensing i?-modes. 



1985). Inflation predicts a nearly scale-invariant spectrum of gravitational waves. 
Their amplitude depends strongly on the energy scale of inflation,'' (power oc Ef 
Rubakov et al 1982; Fabbri & Pollock 1983) and its relationship to the curvature 
fluctuations discriminates between particular models for inflation. Detection of 
gravitational waves in the CMB therefore provides our best hope to study the 
particle physics of inflation. 

Gravitational waves, like scalar fields, obey the Klein-Gordon equation in a flat 
universe and their amplitudes begin oscillating and decaying once the perturba- 
tion crosses the horizon. While this process occurs even before recombination, 
rapid Thomson scattering destroys any quadrupole anisotropy that develops (see 
§3.6). This fact dicates the general structure of the contributions to the power 
spectrum (see Figure 4, left panel): they are enhanced at £ = 2 the present 
quadrupole and sharply suppressed at multipole larger than that of the first peak 
(Abbott & Wise 1984; Starobinsku 1985; Crittenden et al 1993). As is the case 
for the ISW effect, confinement to the low multipoles means that the isolation of 
gravitational waves is severely limited by cosmic variance. 

The signature of gravitational waves in the polarization is more distinct. Be- 
cause gravitational waves cause a quadrupole temperature anisotropy at the end 
of recombination, they also generate a polarization. The quadrupole generated by 
a gravitational wave has its main angular variation transverse to the wavevector 
itself (Hu & White 1997a). The resulting polarization that results has compo- 
nents directed both along or orthogonal to the wavevector and at 45° degree 
angles to it. Gravitational waves therefore generate a nearly equal amount of E 
and B mode polarization when viewed at a distance that is much greater than 

^Ef = V{<p), the potential energy density associated with the scalar field(s) driving inflation. 
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a wavelength of the fluctuation (Kamionkowski et al 1997; Zaldarriaga & Seljak 
1997). The _B-coniponent presents a promising means of measuring the gravita- 
tional waves from inflation and hence the energy scale of inflation (sec Figure 4, 
right panel). Models of inflation correspond to points in the n, Ei plane (Dodel- 
son et al 1997). Therefore, the anticipated constraints will discriminate among 
different models of inflation, probing fundamental physics at scales well beyond 
those accessible in accelerators. 

4.2.4 Gravitational Lensing The gravitational potentials of large-scale 

structure also lens the CMB photons. Since lensing conserves surface bright- 
ness, it only affects anisotropies and hence is second order in perturbation theory 
(Blanchard & Schneider 1987). The photons are deflected according to the an- 
gular gradient of the potential projected along the line of sight with a weighting 
of 2{D^ — D)/(Dt,D). Again the cancellation of parallel modes implies that it is 
mainly the large-scale potentials that are responsible for deflections. Specifically, 
the angular gradient of the projected potential peaks at a multipole £ ^ 60 cor- 
responding to scales of a fc ~ few 10~^ Mpc~^ (Hu 2000b). The deflections are 
therefore coherent below the degree scale. The coherence of the deflection should 
not be confused with its rms value which in the model of Plate 1 has a value of 
a few arcminutes. 

This large coherence and small amplitude ensures that linear theory in the 
potential is sufScient to describe the main effects of lensing. Since lensing is a 
one-to-one mapping of the source and image planes it simply distorts the images 
formed from the acoustic oscillations in accord with the deflection angle. This 
warping naturally also distorts the mapping of physical scales in the acoustic 
peaks to angular scales §3.8 and hence smooths features in the temperature and 
polarization (Seljak 1996a). The smoothing scale is the coherence scale of the 
deflection angle Ai « 60 and is sufficiently wide to alter the acoustic peaks with 
A£ ~ 300. The contributions, shown in Plate 5a are therefore negative (dashed) 
on scales corresponding to the peaks. 

For the polarization, the remapping not only smooths the acoustic power spec- 
trum but actually generates _B-mode polarization (see Plate 1 and Zaldarriaga 
& Seljak 1998). Remapping by the lenses preserves the orientation of the po- 
larization but warps its spatial distribution in a Gaussian random fashion and 
hence does not preserve the symmetry of the original _E-mode. The _B-modes 
from lensing sets a detection threshold for gravitational waves for a finite patch 
of sky (Hu 2001b). 

Gravitational lensing also generates a small amount of power in the anisotropies 
on its own but this is only noticable beyond the damping tail where diffusion has 
destroyed the primary anisotropies (see Plate 5). On these small scales, the 
anisotropy of the CMB is approximately a pure gradient on the sky and the 
inhomogeneous distribution of lenses introduces ripples in the gradient on the 
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scale of the lenses (Seljak & Zaldarriaga 2000). In fact the moving halo effect of 
§4.2.2 can be described as the gravitational lensing of the dipole anisotropy due 
to the peculiar motion of the halo (Birkinshaw & GuU 1983). 

Because the lensed CMB distribution is not linear in the fluctuations, it is 
not completely descril)cd by changes in the power spectrum. Much of the recent 
work in the literature has been devoted to utilizing the non-Gaussianity to isolate 
lensing effects (Bernardeau 1997, 1998; Zaldarriaga & Seljak 1999; Zaldarriaga 

2000) and their cross-correlation with the ISW effect (Goldberg &: Spergel 1999; 
Seljak & Zaldarriaga 1999). In particular, there is a quadratic combination of the 
anisotropy data that optimally reconstructs the projected dark matter potentials 
for use in this cross-correlation (Hu 2001c). The cross correlation is especially 
important in that in a flat universe it is a direct indication of dark energy and 
can be used to study the properties of the dark energy beyond a simple equation 
of state (Hu 2001b). 

4-3 Scattering Secondaries 

From the observations both of the lack of of a Gunn-Peterson trough (Gunn 
& Peterson 1965) in quasar spectra and its preliminary detection (Becker et al 

2001) , wo know that hydrogen was reionized at <^ 6. This is thought to 
occiu' through the ionizing radiation of the flrst generation of massive stars (see 
e.g. Loeb & Barkana 2001 for a review). The consequent recoupling of CMB 
photons to the baryons causes a few percent of them to be rescattered. Linearly, 
rescattering induces three changes to the photon distribution: suppression of 
primordial anisotropy, generation of large angle polarization, and a large angle 
Doppler effect. The latter two are suppressed on small scales by the cancellation 
highlighted in §4.2.1. Non-linear effects can counter this suppression; these are 
the subject of active research and are outlined in §4.3.4. 

4.3.1 Peak Suppression Like scattering before recombination, scattering at 
late times suppresses anisotropies in the distribution that have already formed. 
Reionization therefore suppresses the amplitude of the acoustic peaks by the 
fraction of photons rescattered, approximately the optical depth ~ (see Plate 
5b, dotted line and negative, dashed line, contributions corrcsijonding to |5Aj,|^/'^ 
between the Zri = 7 and Zn = models). Unlike the plasma before recombination, 
the medium is optically thin and so the mean free path and diffusion length of 
the photons is of order the horizon itself. New acoustic oscillations cannot form. 
On scales approaching the horizon at reionization, inhomogeneities have yet to 
be converted into anisotropies (see §3.8) and so large angle fluctuations are not 
suppressed. While these effects are relatively large compared with the expected 
precision of future experiments, they mimic a change in the overall normalization 
of fluctuations except at the lowest, cosmic variance limited, multipoles. 
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4.3.2 Large-Angle Polarization The rescattered radiation becomes po- 
larized since, as discussed in §3.8 temperature inhomogeneities, become anisotro- 
pies by projection, passing through quadrupole anisotropies when the perturba- 
tions are on the horizon scale at any given time. The result is a bump in the power 
spectrum of the _E-polarization on angular scales corresponding to the horizon 
at reionization (see Plate 1). Because of the low optical depth of reionization 
and the finite range of scales that contribute to the quadrupole, the polarization 
contributions are on the order of tenths of //K on scales of ^ ~ few. In a perfect, 
foreground free world, this is not beyond the reach of experiments and can be 
used to isolate the reionization epoch (Hogan et al 1982; Zaldarriaga et al 1997). 
As in the ISW effect, cancellation of contributions along the line of sight guar- 
antees a sharp suppression of coritributious at higher multipolcs in linear theory. 
Spatial modulation of the optical depth due to density and ionization (see §4.3.4) 
does produce higher order polarization but at an entirely negligible level in most 
models (Hu 2000a). 

4.3.3 DoPPLER Effect Naively, velocity fields of order v ~ 10"'^ (sec e.g. 
Strauss & Willick 1995 for a review) and optical depths of a few percent would 
imply a Doppler effect that rivals the acoustic peaks themselves. That this is 
not the case is the joint consequence of the cancellation described in §4.2.1 and 
the fact that the acoustic peaks are not "Doppler peaks" (see §3.8). Since the 
Doppler effect comes from the peculiar velocity along the line of sight, it retains 
no contributions from linear modes with wavevectors perpendicular to the line of 
sight. But as we have seen, these are the only modes that survive cancellation 
(see Plate 3 and Kaiser 1984). Consequently, the Doppler effect from reionization 
is strongly suppressed and is entirely negligible below £ ~ 10^ unless the optical 
depth in the reionization epoch approaches unity (see Plate 5b). 

4.3.4 Modulated Doppler Effects The Doppler effect can survive can- 
cellation if the optical depth has modulations in a direction orthogonal to the 
bulk velocity. This modulation can be the result of either density or ionization 
fluctuations in the gas. Examples of the former include the effect in clusters, and 
linear as well as non-linear large-scale structures. 

Cluster Modulation: The strongly non-linear modulation provided by the 
presence of a galaxy cluster and its associated gas leads to the kinetic Sunyaev- 
Zel'dovich effect. Cluster optical depths on order 10~^ and peculiar velocities 
of 10~^ imply signals in the 10~^ regime in individual arcminute-scale clusters, 
which are of course rare objects. While this signal is reasonably large, it is 
generally dwarfed by the thermal Sunyaev-Zel'dovich effect (see §4.3.5) and has 
yet to be detected with high significance (see Carlstrom et al 2001 and references 
therein). The kinetic Sunyaev-Zel'dovich effect has negligible impact on the power 
spectrum of anisotropies due to the rarity of clusters and can be included as part 
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of the more general density modulation. 

Linear Modulation: At the opposite extreme, linear density fluctuations mod- 
ulate the optical depth and give rise to a Doppler efi'ect as pointed out by Ostriker 
& Vishniac (1986) and calculated by Vishniac (1987) (see also Efstathiou & Bond 
1987). The result is a signal at the fiK level peaking at i ^ few xlO'^ that in- 
creases roughly logarithmically with the reionization redshift (see Plate 5b). 
General Density Modulation: Both the cluster and linear modulations are 
limiting cases of the more general effect of density modulation by the large scale 
structure of the Universe. For the low reionization redshifts currently expected 
(«ri « 6 — 7) most of the effect comes neither from clusters nor the linear regime 
but intermediate scale dark matter halos. An upper limit to the total effect 
can be obtained by assuming the gas traces the dark matter (Hu 2000a) and 
implies signals on the order of ~ few /iK at £ > 10'^ (see Plate 5h). Based on 
simulations, this assumption should hold in the outer proflles of halos (Pearce et al 
2001; Lewis et al 2000) but gas pressure will tend to smooth out the distribution 
in the cores of halos and reduce small scale contributions. In the absence of 
substantial cooling and star formation, these net effects can be modeled under 
the assumption of hydrostatic ocjuilibrium (Komatsu & Scljak 2001) in the halos 
and included in a halo approach to the gas distribution (Cooray 2001). 
Ionization modulation: Finally, optical depth modulation can also come from 
variations in the ionization fraction (Aglianim et al 1996; Gruzinov & Hu 1998; 
Knox et al 1998). Predictions for this effect are the most uncertain as it involves 
both the formation of the first ionizing objects and the subsequent radiative 
transfer of the ionizing radiation (Bruscoli et al 2000; Benson et al 2001). It 
is however unlikely to dominate the density modulated effect except perhaps at 
very high multipoles £ ~ 10^ (crudely estimated, following Gruzinov & Hu 1998, 
in Plate 5b). 

4.3.5 Sunyaev-Zel'dovich Effect Internal motion of the gas in dark mat- 
ter halos also give rise to Doppler shifts in the CMB photons. As in the linear 

Doppler effect, shifts that are first order in the velocity are canceled as photons 
scatter off of electrons moving in different directions. At second order in the ve- 
locity, there is a residual effect. For clusters of galaxies where the temperature of 
the gas can reach Tg ~ lOkeV, the thermal motions are a substantial fraction of 
the speed of light Vrms = (3Te/me)^/^ ~ 0.2. The second order effect represents 
a net transfer of energy between the hot electron gas and the cooler CMB and 
leaves a spectral distortion in the CMB where photons on the Rayleigli- Jeans 
side are transferred to the Wien tail. This effect is called the thermal Sunyaev- 
Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 1972). Because the net effect is 
of order TdmsterTe/me oc UeTg, it is a probe of the gas pressure. Like all CMB 
effects, once imprinted, distortions relative to the rcdsliifting background tem- 
perature remain unaffected by cosmological dimming, so one might hope to find 
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clusters at high redshift using the SZ effect. However, the main effect comes from 
the most massive clusters because of the strong temperature weighting and these 
have formed only recently in the standard cosmological model. 

Great strides have recently been made in observing the SZ effect in individual 
clusters, following pioneering attempts that spanned two decades (Birkinshaw 
1999). The theoretical basis has remained largely unchanged save for small rel- 
ativistic corrections as Te/nie approches unity. Both developemcnts are com- 
prehensively reviewed in (Carlstrom et al 2001). Here we instead consider its 
implications as a source of secondary anisotropies. 

The SZ effect from clusters provides the most substantial contribution to tem- 
perature anisotropies beyond the damping tail. On scales much larger than an 
arcminute where clusters are unresolved, contributions to the power spectrum 
appear as uncorrelated shot noise (Q = const, or At oc i). The additional 
contribution due to the spatial correlation of clusters turns out to be almost neg- 
ligible in comparison due to the rarity of clusters (Komatsu & Kitayama 1999). 
Below this scale, contributions turn over as the clusters become resolved. Though 
there has been much recent progress in simulations (Refregier et al 2000; Seljak 
et al 2001; Springel et al 2001) dynamic range still presents a serious limitation. 

Much recent work has been devoted to semi-analytic modeling following the 
technique of (Cole & Kaiser 1988), where the SZ correlations are described in 
terms of the pressure profiles of clusters, their abundance and their spatial cor- 
relations [now commonly referred to an application of the "halo model" see Ko- 
matsu & Kitayama 1999; Atrio-Barandela & Miicket 1999; Cooray 2001; Komatsu 
& Seljak 2001] . We show the predictions of a simplified version in Plate 5b, where 
the pressure profile is approximated by the dark matter haloprofile and the virial 
temperature of halo. While this treatment is comparatively crude, the inaccu- 
racies that result arc dwarfed by "missing physics"' in both the simulations and 
more sophisticated modelling, e.g. the non-gravitational sources and sinks of en- 
ergy that change the temperature and density profile of the cluster, often modeled 
as a uniform "preheating" of the iutercluster medium (Holder & Carlstrom 2001). 

Although the SZ cfi'cct is expected to dominate the power spectrum of sec- 
ondary anisotropics, it does not necessarily make the other secondaries unmea- 
surable or contaminate the acoustic peaks. Its distinct frequency signature can 
be used to isolate it from other secondaries (see e.g. Cooray et al 2000). Addi- 
tionally, it mainly comes from massive clusters which are intrinsically rare. Hence 
contributions to the power spectrum are non-Gaussian and concentrated in rare, 
spatially localized regions. Removal of regions identified as clusters through X- 
rays and optical surveys or ultimately high resolution CMB maps themselves can 
greatly reduce contributions at large angular scales where they are unresolved 
(Persi et al 1995; Komatsu & Kitayama 1999). 
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4-4 Non-Gaussianity 

As we have seen, most of the secondary anisotropies are not linear in nature and 
hence produce non-Gaussian signatures. Non-Gaussianity in the lensing and SZ 
signals will be important for their isolation. The same is true for contaminants 
such as galactic and extragalactic foregrounds. Finally the lack of an initial non- 
Gaussianity in the fluctuations is a testable prediction of the simplest inflationary 
models (Guth & Pi 1985; Bardeen et al 1983). Consequently, non-Gaussianity 
in the CMB is currently a very active field of research. The primary challenge 
in studies of non-Gaussianity is in choosing the statistic that quantifies it. Non- 
Gaussianity says what the distribution is not, not what it is. The secondary 
challenge is to optimize the statistic against the Gaussian "noise" of the primary 
anisotropies and instrumental or astrophysical systematics. 

Early theoretical work on the bispectrum, the harmonic analogue of the three 
point function addressed its detectability in the presence of the cosmic variance 
of the Gaussian fluctuations (Luo 1994) and showed that the inflationary contri- 
bution is not expected to be detectable in most models (Allen et al 1987; Falk 
et al 1993). The bispectrum is defined by a triplet of multipoles, or configura- 
tion, that defines a triangle in harmonic space. The large cosmic variance in an 
individual configuration is largely offset by the great number of possible triplets. 
Interest was spurred by reports of significant signals in specific bispectrum con- 
figurations in the COBE maps (Ferreira et al 1998) that turned out to be due 
to systematic errors (Banday et al 2000). Recent investigations have focusscd on 
the signatures of secondary anisotropies (Goldberg & Spergel 1999; Cooray & Hu 

2000) . These turn out to be detectable with experiments that have both high 
resolution and angular dynamic range but require the measurement of a wide 
range of configurations of the bispectrum. Data analysis challenges for measur- 
ing the full bispectrum largely remain to be addressed (c.f. Heavens 1998; Spergel 
& Goldberg 1999; Phillips & Kogut 2001). 

The trispectrum, the harmonic analogue of the four point function, also has 
advantages for the study of secondary anisotropies. Its great number of config- 
urations are specified by a quintuplet of multipoles that correspond to the sides 
and diagonal of a quadrilateral in harmonic space (Hu 2001a). The trispectrum is 
important in that it quantifies the covariance of the power spectrum across mul- 
tipoles that is often very strong in non-linear effects, e.g. the SZ effect (Cooray 

2001) . It is also intimately related to the power spectra of quadratic combina- 
tions of the temperature field and has been applied to study gravitational lensing 
effects (Bernardeau 1997; Zaldarriaga 2000; Hu 2001a). 

The bispectrum and trispectrum quantify non-Gaussianity in harmonic space, 
and have clear applications for secondary anisotropies. Tests for non-Gaussianity 
localized in angular space include the Minkowski functionals (including the genus) 
(Winitzki & Kosowsky 1997), the statistics of temperature extrema (Kogut et al 
1996), and wavelet coefficients (Aghanim & Forni 1999). These may be more 
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time-ordered data map(s) bandpowers parameters 

Figure 5: Data pipeline and radical compression. Map are constructed for each frequency 
channel from the data timestreams, combined and cleaned of foreground contamination by 
spatial (represented here by excising the galaxy) and frequency information. Bandpowers are 
extracted from the maps and cosmological parameters from the bandpowers. Each step involves 
a substantial reduction in the number of parameters needed to describe the data, from potentially 
10^° 10 for the Planck sateUite. 



useful for examining foreground contamination and trace amounts of topological 
defects. 

5 DATA ANALYSIS 

The very large CMB data sets that have begun arriving require new, innovative 
tools of analysis. The fundamental tool for analyzing CMB data - the likelihood 
function - has been used since the early days of anisotropy searches (Readhead 
et al 1989; Bond et al 1991; Dodelson & Jubas 1993). Brute force likehhood 
analyses (Tegmark & Bunn 1995) were performed even on the relatively large 
COBE data set, with six thousand pixels in its map. Present data sets are a 
factor of ten larger, and this factor will soon get larger by yet another factor of 
a hundred. The brute force approach, the time for which scales as the number of 
pixels cubed, no longer suffices. 

In response, analysts have devised a host of techniques that move beyond the 
early brute force approach. The simplicity of CMB physics ~ due to linearity 
- is mirrored in analysis by the apparent Gaussianity of both the signal and 
many sources of noise. In the Gaussian limit, optimal statistics are easy to iden- 
tify. These compress the data so that all of the information is retained, but the 
subsequent analysis - because of the compression - becomes tractable. 

The Gaussianity of the CMB is not shared by other cosmological systems since 
gravitational non-linearities turn an initially Gaussian distribution into a non- 
Gaussian one. Nontheless, many of the techniques devised to study the CMB have 
been proposed for studying: the 3D galaxy distribution (Tegmark et al 1998), the 
2D galaxy distribution (Efstathiou & Moody 2001; Huterer et al 2001) the Lyman 
alpha forest (Hui et al 2001), the shear field from weak lensing (Hu & White 
2001), among others. Indeed, these techniques are now indispensible, powerful 
tools for all cosniologists, and we would be remiss not to at least outline them in 
a disussion of the CMB, the context in which many of them were developed. 
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Figure 5 summarizes the path from the data analyis starting point, a timestream 
of data points, to the end, the determination of cosmological parameters. Pre- 
ceding this starting point comes the calibration and the removal of systematic 
errors from the raw data, but being experiment specific, we do not attempt to 
cover such issues here.'* Each step radically compresses the data by reducing 
the number of parameters used to describe it. Although this data pipeline and 
our discussion below are focused on temperature anisotropies, similar steps have 
been elucidated for polarization (Bunn 2001; Tegmark & de Oliveira-Costa 2001; 
Lewis et al 2001). 



5. 1 Mapmaking 



An experiment can be characterized by the data dt taken at many different times; 
a pointing matrix Pti, relating the data timestream to the underlying signal at 
pixelized positions indexed by i, and a noise matrix C^^tf characterizing the 
covariance of the noise in the timestream. A model for the data then is dt = 
Pti'S>i + rit (with implicit sum over the repeating index i); it is the sum of signal 
plus noise. Here nt is drawn from a distribution (often Gaussian) with mean 
zero and covariance {ritn-ti) = C^^tt'- In its simplest form the pointing matrix P 
contains rows - which corresponds to a particular time ~ with all zeroes in it 
except for one column with a one (see Figure 5). That column corresponds to 
the particular pixel observed at the time of interest. Typically, a pixel will be 
scanned many times during an experiment, so a given column will have many 
ones in it, corresponding to the many times the pixel has been observed. 

Given this model, a well-posed question is: what is the optimal estimator for 
the signal 0i? i.e. what is the best way to construct a map? The answer stems 
from the likelihood function £, defined as the probability of getting the data given 
the theory £ = P[data| theory]. In this case, the theory is the set of parameters 



£e{dt 



(27r)^*/2VditC: 



; exp 



(26) 



That is, the noise, the difference between the data and the modulated signal, is 
assumed to be Gaussian with covariance Cd- 

There are two important theorems useful in the construction of a map and 
more generally in each step of the data pipeline (Tegmark et al 1997). The first 
is Bayes' Theorem. In this context, it says that P[0i|dt], the probability that 
the temperatures are equal to 0.; given the data, is proportional to the likelihood 
function times a prior P{Qi). Thus, with a uniform prior. 



P[0,;|dt] (X P[dtm = Ceidt 



Aside from COBE, experiments to date have had a sizable calibration error 
must be factored into the interpretation of Plate 1. 



(27) 

' 5-10%) which 
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with the nomiahzation constant determined by requiring the integral of the prob- 
ability over all Qi to be equal to one. The probability on the left is the one of 
interest. The most likely values of Qi therefore are those which maximize the 
likelihood function. Since the log of the likelihood function in question, Equation 
(26), is quadratic in the parameters 0^, it is straightforward to find this maxi- 
mum point. Differentiating the argument of the exponential with respect to Qi 
and setting to zero leads immediately to the estimator 

©i = CN,ijPjtCju,dt' , (28) 

where Cat = (P*^C^^P)~^. As the notation suggests, the mean of the estimator 
is equal to the actual 0; and the variance is ecjual to Ca?. 

The second theorem states that this maximum likelihood estimator is also the 
minimum variance estimator. The Cramer-Rao inequality says no estimator can 
measure the Gj with errors smaller than the diagonal elements of F~^, where the 
Fisher matrix is defined as 

--(-SI)- 

Inspection of Equation (26) shows that, in this case the Fisher matrix is precisely 
equal to C^'^ . Therefore, the Cramer-Rao theorem implies that the estimator of 
Equation (28) is optimal: it has the smallest possible variance (Tegmark 1997a). 
No information is lost if the map is used in subsequent analysis instead of the 
timestream data, but huge factors of compression have been gained. For example, 
in the recent Boomerang experiment (Netterfield et al 2001), the timestream 
contained 2 x 10* numbers, while the map had only 57, 000 pixels. The map 
resulted in compression by a factor of 3500. 

There are numerous complications that must be dealt with in realistic appli- 
cations of Equation (28). Perhaps the most difficult is estimation of C,;, the 
timestream noise covariance. This typically must be done from the data itself 
(Ferreira & Jaffe 2000; Stompor et al 2001). Even if were known perfectly, 
evaluation of the map involves inverting C^, a process which scales as the number 
of raw data points cubed. For both of these problems, the assumed stationarity of 
Cd,tt' (it depends only on f — t') is of considerable utility. Iterative techniques to 
approximate matrix inversion can also assist in this process (Wright et al 1996). 
Another issue which has received much attention is the choice of pixelization. 
The community has converged on the Healpix pixelization scheme^, now freely 
available. 

Perhaps the most dangerous complication arises from astropliysical foregrounds, 
both within and from outside the Galaxy, the main ones being synchrotron, 
bremmsstrahlung, dust and point source emission. All of the main foregrounds 

^http ; / /wwu . eso . org/science/healpix/ 
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have different spectral shapes than the blackbody shape of the CMB. Modern 
experiments typically observe at several different frequencies, so a well-posed 
question is: how can we best extract the CMB signal from the different fre- 
quency channels (Bouchet & Gispert 1999)? The blackbody shape of the CMB 
relates the signal in all the channels, leaving one free parameter. Similarly, if 
the foreground shapes are known, each foreground comes with just one free pa- 
rameter per pixel. A likelihood function for the data can again be written down 
and the best estimator for the CMB amplitude determined analytically. While in 
the absence of foregrounds, one would extract the CMB signal by weighting the 
frequency channels according to inverse noise, when foregrounds arc present, the 
optimal combination of different frequency maps is a more clever weighting that 
subtracts out the foreground contribution (Dodelson 1997). One can do better if 
the pixel-to-pixel correlations of the foregrounds can also be modeled from power 
spectra (Tegmark & Efstathiou 1996) or templates derived from external data. 

This picture is complicated somewhat because the foreground shapes are not 
precisely known, varying across the sky, e.g. from a spatially varying dust tem- 
perature. This too can be modelled in the covariance and addressed in the likeli- 
hood (Tegmark 1998; White 1998). The resulting cleaned CMB map is obviously 
noisier than if foregrounds were not around, but the multiple channels keep the 
degradation managable. For example, the errors on some cosmological parame- 
ters coming from Planck may degrade by almost a factor of ten as compared with 
the no-foreground case. However, many errors will not degrade at all, and even 
the degraded parameters will still be determined with unprecedented precision 
(Knox 1999; Prunet et al 2000; Tegmark et al 2000). 

Many foregrounds tend to be highly non-Gaussian and in particular well- 
localized in particular regions of the map. These pixels can be removed from 
the map as was done for the region around the galactic disk for COBE. This 
technique can also be highly effective against point sources. Indeed, even if there 
is only one frequency channel, external foreground templates set the form of 
the additional contributions to Cat. which, when properly included, immunize 
the remaining operations in the data pipeline to such contaminants (Bond et al 
1998). The same technique can be used with templates of residual systematics 
or constraints imposed on the data, from e.g. the removal of a dipole. 

5.2 Bandpower Estimation 

Figure 5 indicates that the next step in the compression process is extracting 
bandpowers from the map. What is a bandpower and how can it be extracted from 
the map? To answer these questions, we must construct a new likelihood function, 
one in which the estimated Qi are the data. No theory predicts an individual 
Qi, but all predict the distribution from which the individual temperatures are 
drawn. For example, if the theory predicts Gaussian fluctuations, then is 
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distributed as a Gaussian with mean zero and covariance equal to tlie sum of tlie 
noise covariance matrix Cn and the covariance due to the finite sample of the 
cosmic signal Cg. Inverting Equation (1) and using Equation (2) for the ensemble 
average leads to 

Cs,ij = {Oiej) = A^^^Wi^ij , (30) 
e 

where A^,^ depends on the theoretical parameters through Ce (see Equation 
(3)). Here We., the window function, is proportional to the Legendre polynomial 
Pe{hi ■ Uj) and a beam and pixel smearing factor 6|. For example, a Gaussian 
beam of width a, dictates that the observed map is actually a smoothed picture 
of true signal, insensitive to structure on scales smaller than a. If the pixel scale 
is much smaller than the beam scale, 6| cx , Techniques for handling 

asymmetric beams have also recently been developed (Wu ct al 2001; Wandelt 
& Gorski 2001; Souradeep & Ratra 2001). Using bandpowers corresponds to 
assuming that A^^ is constant over a finite range, or band, of £, equal to Ba for 

— '5^o/2 < i < ia + 5£a/2. Plate 1 gives a sense of the width and number of 
bands AT;, probed by existing experiments. 

For Gaussian theories, then, the likelihood function is 

- (2.)^wVdetCe ' ^''^ 

where Ce = C5 + and Np is the number of pixels in the map. As before, 
£b is Gaussian in the anisotropics Bj, but in this case Oi are not the parameters 
to be determined; the theoretical parameters are the Ba, upon which the covari- 
ance matrix depends. Therefore, the likelihood function is not Gaussian in the 
parameters, and there is no simple, analytic way to find the point in parameter 
space (which is multi-dimensional depending on the number of bands being fit) at 
which £b is a maximum. An alternative is to evaluate £b numerically at many 
points in a grid in parameter space. The maximum of the Cb on this grid then 
determines the best fit values of the parameters. Gonfidence levels on say Bi can 
be determined by finding the region within which J^dBi[n^22ldf]CB = 0.95, 
say, for 95% limits. 

This possibility is no longer viable due to the sheer volume of data. Consider 
the Boomerang experiment with Np = 57, 000. A single evaluation of £b involves 
computation of the inverse and determinant of the Np x Np matrix Ce, both 
of which scale as N^. While this single evaluation might be possible with a 
powerful computer, a single evaluation does not suffice. The parameter space 
consists of Ni, = 19 bandpowers equally spaced from la = 100 up to la = 1000. 
A blindly placed grid on this space would require at least ten evaluations in each 
dimension, so the time required to adequately evaluate the bandpowers would 
scale as 10"'^'' A''^. No computer can do this. The situation is rapidly getting worse 
(better) since Planck will have of order 10*^ pixels and be sensitive to of order a 
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10-* bands. 

It is clear that a "smart'" sampling of the likelihood in parameter space is nec- 
essary. The numerical problem, searching for the local maximum of a function, 
is well-posed, and a number of search algorithms might be used. Cb tends to 
be sufficiently structureless that these techniques suffice. Bond et al (1998) pro- 
posed the Newton-Raphson method which has become widely used. One expands 
the derivative of the log of the likelihood function - which vanishes at the true 
maximum of £b - around a trial point in parameter space, Ba'\ Keeping terms 
second order in Ba — si^' leads to 

4 = Bf (32) 

where the curvature matrix ^6 is the second derivative of — In £b with respect 
to Ba and Note the subtle distinction between the curvature matrix and 
the Fisher matrix in Equation (29), F = (F). In general, the curvature matrix 
depends on the data, on the 0;. In practice, though, analysts typically use the 
inverse of the Fisher matrix in Equation (32). In that case, the estimator becomes 

Ba = + (e,Ce,i.^Ce,Le, - C^l,^) , (33) 

quadratic in the data 0j. The Fisher matrix is equal to 

In the spirit of the Newton-Raphson method. Equation (33) is used itcrativcly 
but often converges after just a handful of iterations. The usual approximation 
is then to take the covariance between the bands as the inverse of the Fisher 
matrix evaluated at the convergent point Cb = F^^. Indeed, Tegmark (1997b) 
derived the identical estimator by considering all unbiased quadratic estimators, 
and identifying this one as the one with the smallest variance. 

Although the estimator in Equation (33) represents a 10^'" improvement 
over brute force coverage of the parameter space - converging in just several 
iterations - it still requires operations which scale as Np. One means of speeding 
up the calculations is to transform the data from the pixel basis to the so-called 
signal-to-noise basis, based on an initial guess as to the signal, and throwing out 
those modes which have low signal-to-noise (Bond 1995; Bunn & Sugiyama 1995). 
The drawback is that this procedure still requires at least one A^^ operation and 
potentially many as the guess at the signal improves by iteration. Methods to 
truly avoid this prohibitive A''^ scaling (Oh et al 1999; Wandelt Sz Hansen 2001) 
have been devised for experiments with particular scan strategies, but the general 
problem remains open. A potentially promising approach involves extracting the 
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real space correlation functions as an intermediate step between the map and 
the bandpowers (Szapudi et al 2001). Another involves consistently analyzing 
coarsely pixelized maps with finely pixelized sub-maps (Dore et al 2001). 



5.3 Cosmological Parameter Estimation 

Tlic huge advantage of bandpowers is that they represent the natural meeting 
ground of theory and experiment. The above two sections outline some of the 
steps involved in extracting them from the observations. Once they are extracted, 
any theory can be compared with the observations without knowledge of exper- 
imental details. The simplest way to estimate the cosmological parameters in a 
set Cj is to approximate the likelihood as 



1 



(27r)^<^/VdetCB 



exp 



-{Ba-Ba)CgUBt-Bt) 



(35) 



and evaluate it at many points in parameter space (the bandpowers depend on 
the cosmological parameters). Since the number of cosmological parameters in 

the working model is Nc ^ 10 this represents a final radical compression of 
information in the original timestrcam which recall has up to Nt ~ 10^" data 
points. 

In the approximation that the band power covariance C_b is independent of 
the parameters c, maximizing the likelihood is the same as minimizing x^- This 
has been done by dozens of groups over the last few years especially since the 
release of CMBFAST (Seljak & Zaldarriaga 1996), which allows fast computation 
of theoretical spectra. Even after all the compression summarized in Figure 5, 
these analyses are still computationally cumbersome due to the large numbers 
of parameters varied. Various methods of speeding up spectra computation have 
been proposed (Tegmark & Zaldarriaga 2000), based on the imderstanding of the 
physics of peaks outlined in §3, and Monte Carlo explorations of the likelihood 
function (Christenscii ct al 2001). 

Again the inverse Fisher matrix gives a quick and dirty estimate of the errors. 
Here the analogue of Equation (29) for the cosmological parameters becomes 



dCi dCj 



(36) 



In fact, this estimate has been widely used to forecast the optimal errors on cos- 
mological parameters given a proposed experiment and a band covariance matrix 
Cb which includes diagonal sample and instrumental noise variance. The reader 
should be aware that no experiment to date has even come close to achieving the 
precision implied by such a forecast! 

As wc enter the age of precision cosmology, a number of caveats will become in- 
creasingly important. No theoretical spectra are truly fiat in a given band, so the 
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question of how to weight a theoretical spectrum to obtain Ba can be important. 
In principle, one must convolve the theoretical spectra with window functions 
(Knox 1999) distinct from those in Equation (30) to produce Ba- Among re- 
cent experiments, DASI (Pryke et al 2001) among others have provided these 
functions. Another complication arises since the true likelihood function for Ba 
is not Gaussian, i.e. not of the form in Equation (35). The true distribution is 
skewed: the cosmic variance of Equation (4) leads to larger errors for an upward 
fiuctuation than for a downward fiuctuation. The true distribution is closer to 
log-normal (Bond et al 2000), and several groups have already accounted for this 
in their parameter extractions. 

6 DISCUSSION 

Measurements of the acoustic peaks in the CMB temperature spectrum have 
already shown that the Universe is nearly spatially fiat and began with a nearly 
scale-invariant spectrum of curvature fluctuations, consistent with the simplest 
of inflationary models. In a remarkable confirmation of a longstanding prediction 
of Big Bang Nucleosynthesis, the CMB measurements have now verified that 
baryons account for about four percent of the critical density. Further, they 
suggest that the matter density is some ten times higher than this, implying the 
existence of non-baryonic dark matter and dark energy. 

Future measurements of the morphology of the peaks in the temperature and 
polarization should determine the baryonic and dark matter content of the Uni- 
verse with exquisite precision. Beyond the peaks, gravitational wave imprint on 
the polarization, the gravitational lensing of the CMB, and gravitational and 
scattering secondary anisotropies hold the promise of understanding the physics 
of inflation and the impact of dark energy on structure formation. 

The once and future success of the CMB anisotropy enterprise rests on three 
equally important pillars: advances in experimental technique, precision in the- 
ory, and development of data analysis techniques. The remarkable progress in 
the field over the last decade owes much to the efforts of researchers in all three 
disciplines. That much more effort will be required to fulfill the bright promise 
of the CMB suggests that the field will remain active and productive for years to 
come. 
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